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Abstract 
The age-old idea of " like dissolves like" requires a notion of " likeness" that is hard to quantify for polymers. 
We revisit the concepts of pure component cohesive energy density c^ and mutual cohesive energy density 
Ci2 so that they can be extended to polymers. We recognize the inherent limitations of Ci2 due to its very 
Q ' definition, which is based on the assumption of no volume of mixing (true for incompressible systems), one 
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of the assumptions in the random mixing approximation (RMA); no such limitations are present in the 

identification of c^. We point out that the other severe restriction on ci2 is the use of pure components 

, ■ in its definition because of which ci2 is not merely controlled by mutual interactions. Another quantity 

Ch ' cf^® as a measure of mutual cohesive energy density that does not suffer from the above limitations of 

O I Ci2 is introduced. It reduces to Ci2 in the RMA limit. We are able to express cff^^ in terms of ci2 and 

pure component c^'s. We also revisit the concept of the internal pressure and its relationship with the 

conventional and the newly defined cohesive energy densities. In order to investigate volume of mixing 

f^ ' effects, we introduce two different mixing processes in which volume of mixing remains zero. We then carry 

- ' out a comprehensive reanalysis of various quantities using a recently developed recursive lattice theory that 



o 

^D i theory such as the Flory-Huggins theory for polymers. In the RMA limit, our recursive theory reduces to 



has been tested earlier and has been found to be more accurate than the conventional regular solution 



C^ • the Flory-Huggins theory or its extension for a compressible blend. Thus, it supersedes the Flory-Huggins 

theory. Consequently, the conclusions based on our theory are more reliable and should prove useful. 
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I. INTRODUCTION 

A. Pure Component and Athermal Reference State 

Understanding the solubility of a polymer in a solvent is a technologically important problem. 
It is well documented that "like dissolves like" but it is almost impossible to quantify the notion of 
"likeness" of materials. The understanding of solubility requires a basic understanding of "likeness" 
that is lacking at present. Solubility parameters in all their incarnations are attempts to quantify 
this simple notion. There are several ways one can proceed to understand solubility by considering 
various thermodynamic quantities, not all of which are equivalent. However, almost all these 
current approaches are based on our thermodynamic understanding of mixtures at the level of the 
regular solution theory (RST) l|, |2,y|- 

In practice, one usually introduces the Hildebrand solubility parameter 



for a pure component in terms of c^, the latter being the cohesive energy density of the pure (P) 
component and is normally reported at its boiling temperature. In this work, we will take this 
definition to be operational at any temperature, thus treating 5 as a thermodynamic quantity, and 
not just a parameter. Then its value at the boiling point will be the customarily quoted solubility 
parameter. The cohesive energy is related to the interaction energy £[^t obtained by subtracting 
the energy i?ath of the hypothetical athermal state of the pure component, the state in which the 
self- interactions (both interparticle and intraparticle) are absent, from the total energy £: 

"^int = £ — f ath- (1) 

The athermal state is usually taken to be at the same temperature T and the volume V as the 
system itself. In almost all cases of interest, iSath is nothing but the kinetic energy and depends 
only on T but not on V. Thus, iSmt depends directly on the strength of the self-interaction, the only 
interaction present in a pure component and vanishes as self-interactions disappear, since £ — > i?ath 
as this happens. Thus, iSint can be used to estimate the strength of the self-interaction. One can 
also take the hypothetical state to be at the same T and the pressure P. In this case, there would 
in principle be a difference in the volume V of the pure component and T4th of the hypothetical 
state, but this diff'erence will not change iSint in ©• The hypothetical state is approximated in 
practice by the vapor phase at the boiling point in which the particles are assumed to be so far 
apart that their mutual interactions can be neglected. However, to be precise, it should be the 



vapor phase at zero pressure so that the volume is infinitely large to make the particle separation 
also infinitely large to ensure that they are non-interacting. This causes problems for polymers 
[J]. Our choice of the athermal state in ^ to define £^int overcomes these problems altogether. By 
definition, 

cP = -£i^t/V (2) 

is the negative of the interaction energy density (per unit volume V) of the system at a given 
temperature T. (At the boiling point, V is taken to be the volume of the liquid). The negative sign 
is to ensure c^ > since usually S-mt is negative to give cohesion. Because of its dimensions, c^ is 
also known as the cohesive pressure. In this form, c^ is a thermodynamic quantity and represents 
the thermodynamically averaged (potential) energy per unit volume of the pure component. Thus, 
c can be calculated even for macromolecules like polymers, which is of central interest in this 
work or for polar solvents. It is a macroscopic, i.e. a thermodynamic quantity characterizing 
microscopic interparticle interactions in a pure component. This is important as it is well known 
that 5 cannot be measured directly as most polymers cannot be vaporized without decomposing 
[5|. Thus, theoretical means are needed to evaluate 6, which is our goal in this work. 

It should be noted from the above definition that c^ contains the contributions from all kinds of 
forces (van der Waals, dipolar, and hydrogen bonding forces) in the system J6|. In this work, we are 
only interested in the weak van der Waals interactions for simplicity, even though the investigation 
can be easily extended to include other interactions. It should also be noted that the definition (jSJ 
does not suffer from any inherent approximation, and can be used to calculate c in any theory 
or by experimental measurements. As we will see, this is not true of the mutual cohesive density 
definition, which is introduced below. 

B. Mixture and Self-interacting Reference State 

As it stands, the pure component quantity c^ is oblivious to what may happen in a mixture 
formed with another component. Despite this, c^ or 6 is customarily interpreted as a rough measure 
of a solvent's ability to dissolve a solute ("like dissolves like"). This interpretation of the solubility 
parameter is supposed to be reliable only for non-polar solvents formed of small molecules, and one 
usually refrains from using it for polar solvents such as esters, ketones, alcohol, etc. Our interest 
is to investigate this quantity for macromolecules here, and its significance for the solubility in a 
mixture. According to the famous Scatchard-Hildebrand relation 3, the energy of mixing AEm 



per unit volume for a binary mixture of two components 1 and 2 must always be non-negative since 
it is given by 

AEm = (Si - 62f^Pl(p2, {^1 + f2 = 1), (3) 

where (/Jj are the volume fractions of the two components i, and 5i and 82 are their solubility 
parameters. It is implicitly assumed here that the volume of mixing I^Vm = 0. Thus, this expression 
does not contain the contribution from a non-zero volume of mixing. We will be interested in 
investigating this additional contribution in this work. Later, we will discover that Q can only be 
justified [see (|1T|) below], if we take 

<^j = 4>rai/(t>ni-, (</'m = 0ml + 0m2) (4) 

as representing the monomer density ratios, or monomer fractions; see below for precise definition. 
Only in the RST can we identify (pi with the volume fraction of the ith component. The significance 
of (121) is that the behavior of the mixture is completely determined by the pure component properties 
and provides a justification for "like dissolves like". This must be a gross approximation even for 
„„n-p„l« systems and ca„„o. be t™e in genera, since .he energy of mixing can be negative in 
many cases, as shown elsewhere |9|, and as we will also see here later; see, for example. Fig. 1161 
What we discover is that lS.Ey[ can be negative even if AT4i = 0. Thus, zero volume of mixing is 
not sufficient for the validity of (jH)). 

For the mixture, we need to introduce a thermodynamic or macroscopic quantity ci2 character- 
izing the mutual interaction between the two components; this quantity should be determined by 
the microscopic interaction between the two components. Therefore, the introduction of ci2 will, in 
principle, require comparing the mixture with a hypothetical mixture in which the two components 
have no mutual interactions similar to the way the pure component (having self-interaction) was 
compared with the athermal state above (in which the self-interaction was absent) for the intro- 
duction of c^. The hypothetical state of the mixture should not be confused with the athermal 
state of the mixture; the latter will require even the self interaction of the two components to be 
absent. The new hypothetical state will play an important role in our investigation and will be 
called self-interacting reference state (SRS). The mutual interaction energy in the binary mixture 
is obtained by subtracting the energy of SRS from that of the mixture: 

4f ^£-£sKs; (5) 

compare with (^. Just as before, SRS can be taken at the same T and V, or T and P as the 
mixture. This allows us to separate out the two contributions, one due to the presence of mutual 



interactions at the same volume Vsrs = V of SRS, and the second contribution due to merely a 
change in the volume from Vsrs to V; this was not the case with Si^t above. Each contribution of 
£^^^ vanishes as mutual interactions disappear, since £ -^ <?SRS) and Vsrs -^ V as this happens. 
If £'int is used to introduce a mutual cohesive energy density (to be denoted later by cf^^^), then 
such a density would most certainly vanish with vanishing mutual interaction strength. However, 
this is not the conventional approach adopted in the literature when introducing ci2- Rather, one 
considers the energy of mixing. We will compare the two approaches in this work. Whether the 
conventionally defined C12 vanishes with the mutual interactions remains to be verified. In addition, 
whether it is related to the pure component self interaction cohesive energy densities Cii and C22 in 
a trivial fashion such as H1U() . see below, needs to be investigated. As we will see, this will require 
further assumptions even within RST to which we now turn. 

C. Regular Solution Theory (RST) 

The customary approach to introduce C12 is to follow the classical approach developed by van 



Laar, and Hildebrand 
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which is based on RST, a theory that can be developed consistently 
only on a lattice. The theory describes an incompressible lattice model or a compressible model in 
which the volume of mixing AT^vi is zero. The lattice model is introduced as follows. One usually 
takes a homogeneous lattice of coordination number q and containing N sites. The volume of the 
lattice is Nvq where vq is the lattice cell volume. We place on this lattice the polymer (component 
1) and the solvent (component 2) molecules in such a way that the connectivity of the molecules 
are kept intact. It should be stressed that the solvent molecules themselves may be polymeric in 
nature. In addition, the excluded volume interactions are enforced by requiring that no site can 
be occupied by more than one monomer at a time. The monomer densities of the two components 
4'rai {i = 1) 2) are the densities of sites occupied by the ith component. Two monomers belonging 
to components i and j, respectively, interact with an energy Cij = Cji only when they are nearest- 
neighbor. (This nearest-neighbor restriction can be easily relaxed, but we will not do that here for 
simplicity.) The interaction between the polymer and the solvent is then characterized by a single 
excess energy parameter e = e\2 defined in general by the combination 

Eij =eij-(l/2)(eii-hejj). (6) 



The origin of this combination is the fixed lattice connectivity as shown elsewhere J7|. For an 
incompressible binary mixture for which 

(/•ml + 0m2 = 1, 

the excess energy e is sufficient for a complete thermodynamic description on a lattice [3|- 

On the other hand, a compressible lattice model of the mixture, which requires introducing voids 
as a separate component (i = 0), will usually require two additional excess energy parameters eoij 
and £02 |3l- I^ the following, we will implicitly assume that a void occupies a site of the lattice and 
has it occupies a volume vq. In our picture, a pure component is a pseudo-binary mixture (i = 0, 1 
or i = 0, 2), while a compressible binary mixture is a pseudo-tertiary mixture {i = 0, 1, 2). 
Since voids do not interact with themselves or with any other component, we must set 

eoo = eoi = 0, 

so that the corresponding excess energy 

eoi = -(l/2)eii, (7) 

see ©, and is normally positive since en is usually negative. 

Two of the three conditions for RST to be operative are 

(i) Isometric Mixing: AVm = 0, and 

(ii) Symmetric Mixture: The two components must be of the same size and have the same 
interaction (en = 622)- This is a restatement of "like dissolves like." 

The condition Al^i = for isometric mixing is always satisfied in an incompressible system. 
For a compressible system, AVm need not be zero, and one can have either an isometric or a 
non- isometric mixing depending on the conditions. 

1. Random Mixing Approximation (RMA) 

The third condition for RST to be valid is the fulfillment of the 

iii) RMA Limit: The interaction energy Sij be extremely weak {eij — > 0), and the coordination 
number q of the lattice extremely large {q -^ 00) simultaneously so that the product 

qsij remains fixed and finite, as Sij — > 0, and q — > 00. 



Indeed, if one introduces the diniensionless Flory-Huggins chi parameter Xij = PQ^ijy where 
P = l/k^T, T being the temperature in the Kelvin scale, then one can also think of keeping Xij 
fixed and finite, instead of qSij, under the simultaneous limit 



Xij = (^Q^ij fixed and finite, as Psij -^ 0, and q ^ oo. 



(8) 



It is quite useful to think of RST in terms of these limits as both Psij and q are dimensionless. The 
above simultaneous limit gives rise to what is usually known as the random mixing approximation 
(RMA), and has been discussed in detail in a recent review article J9| in the context of polymer 
mixtures. For an incompressible system, we need to keep only a single chi parameter x = Pq^ 
fixed and finite in the limit. For a compressible system, one must also simultaneously keep the two 
additional chi parameters related to Sqi and £02 , and pPvo fixed and finite, where P is the pressure 



The RMA limit can be applied even to a pure component (for which the first two conditions 
are meaningless). It can also be applied when mixing is not isometric [Hl or when the mixture is 
not symmetric. Therefore, RST is equivalent to the isometric RMA. In the unusual case X = Oj 
{q — > cxd), the resulting isometric theory is known as the ideal solution theory, which will not be 
considered here as we are interested in the case x 7^ 0. 



2. London- Berthelot and Scatchard-Hildebrand Conjectures 



The energy of mixing AEm per unit volume is the central 
and can be used to introduce an effective "energetic" chi 
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XeS = PAEMVo/<Pml4>m2, 



uantity for solubility considerations, 
as follows: 

(9) 



which is a measure of the Flory-Huggins x(= X12) parameter or the excess energy e{= £12); the 
latter is directly related to the mutual interaction energy ei2, see ©, which explains the usefulness 
of A£^M for solubility considerations. One of the important consequences of the application of RST 
is the Scatchard-Hildebrand conjecture |3|, according to which the energy of mixing is given by a 
non- negative form Q. This is known to be violated; see [9|, and below. One of the reasons for its 
failure could be the much abused London-Berthelot conjecture 



C12 = \ cf-iC^o = (5152, 



(10) 



relating the mutual cohesive energy density ci2 with the pure component cohesive energy densities 
c^i, C22, and used in the derivation of @. In contrast, the London conjecture 



ei2 



V(-eii)(-e22) (11) 



deals directly with the microscopic interaction energies, and is expected to hold for non-polar 
systems; see also 1,11] for an illuminating derivation and discussion. In the isometric RMA limit, 
the two conjectures become the same in that (|11|) implies (jiup . In general, they are two independent 
conjectures. As we will demonstrate here, (|11() does not imply (|1U() once we go beyond the RMA 
limit. We will also see that the non- negativity of the form ^ is violated even for isometric mixing. 
Association of one component, such as through hydrogen bonding, usually makes 6^2 < 
|eii| 1 622 1 • On the other hand, complexing results in the opposite inequality e^2 > kill 1 622! • It 
is most likely that the binary interaction between monomers of two distinct species will deviate 
from the London conjecture (jllj) to some degree. Thus, some restrictions have to be put on the 
possible relationship between these energy parameters for our numerical calculations. We have 
decided to consider only those parameters that satisfy the London conjecture (jllI) in this work for 
physical systems. 

3. Deviation from London- Berthelot Conjecture 

The deviation from the London-Berthelot conjecture pUj) is usually expressed in terms of a 
binary quantity I12 defined via 



Cl2 = ^f^2{^-h2), (12) 

and it is usually believed that the magnitude of /12 is very small: 

I/12I « 1. (13) 

It is possible that I12 vanishes at isolated points, so that the London-Berthelot conjecture becomes 
satisfied. This does not mean that the system obeys RST there. As we will demonstrate here, 
we usually obtain a non-zero /12, thus implying a failure of the London-Berthelot conjecture (|lfl|>. 
even if the London conjecture (|11() is taken to be operative. Another root for the failure of H10|) 
under this assumption could be non-isometric mixing. A third cause for the failure could be the 
corrections to the RMA limit, since a real mixture is not going to truly follow RST. To separate 
the effects of the three causes, we will pay particular attention to I12 in this work. We will assume 



the London conjecture Ullj) to be valid, and consider the case of isometric mixing. We will then 
evaluate I12 using a theory that goes beyond RST. A non-zero I12 in this case will then be due to 
the corrections to the RMA limiting behavior or RST. 

D. Internal Pressure 

Hildebrand [j] has also argued that the solubility of a given solute in different solvents is 
determined by relative magnitudes of internal pressures, at least for nonpolar fluids. Thus, we will 
also investigate the internal pressure, which is defined properly by 

Pm = -Pint = -{P-P.th), (14) 

where Pint is the contribution to the pressure P due to interactions in the system, and is obtained 
by subtracting Path from P, where Path is the pressure of the hypothetical athermal state. The 
volume V of the system, which has the pressure P, may or may not be equal to the volume of 
the hypothetical athermal state. This means that Pin will have different values depending on the 
choice of athermal state volume. In this work, we will only consider the athermal state whose 
volume is equal to the volume V of the system; its pressure then has the value Path; which is 
not equal to the pressure P of the system. As the interactions appear in the athermal system at 
constant volume, its pressure Path will reduces due to attractive interactions. This will give rise 
to Pint, which is negative so that Pin will be a positive quantity for attractive interactions. For 
repulsive interactions, which we will not consider here. Pin will turn out to be a negative quantity. 
In either case, as we will show here, Pin should be distinguished from {d£/dV)T, with which it is 
usually equated. Their equality holds only in a very special case as we will see here. 

The availability of PyT-data makes it convenient to obtain {d£/dV)T- Therefore, it is not 
surprising to find it common to equate it with Pin- In a careful analysis, Sauer and Dee have 
shown a close relationship between {d£/dV)T and c |l2l ]. Here, we will investigate the relationship 
between Pin and c^. 

E. Internal Volume 

One can also think of an alternate scenario in which the pressure of the athermal state is kept 
constant as interactions appear in the athermal system. This pressure can be taken to be either P 
or Path- Let us keep its pressure to be P, and so that its volume is T4th- In this case, the volume 
of the system V will be smaller (greater) than Fath of the athermal state because of the attractive 



(repulsive) interactions, so that one can also use the negative of the change in the volume V — T4th 
as a measure of the attractive (repulsive) interactions. This allows us to introduce the following 
volume, to be called internal volume 

Vm = -Mnt = -{V - Fath) (15) 

as a measure of cohesiveness. 

The two internal quantities are mutually related and one needs to study only one of them. 

F. Beyond Regular Solution Theory: Solubility and Effective Chi 

In polymer solutions or blends of two components 1 and 2, non-isometry is a rule rather than 
exception due to asymmetry in the degrees of polymerization and in the pure component interac- 
tions. Thus, the regular solution theory is most certainly inoperative. An extension of RST, to be 
described as the non- isometric regular solution theory, allows for non-isometric mixing (AVm t^ 0), 
and its successes and limitations have been considered elsewhere |l3 . where it is simply called the 
regular solution theory]. It is the extended theory that will be relevant for polymers. 

Let us suppress the pure component index i in the following. Crudely speaking, c^ for a 
pure component is supposed to be related to the pure component interaction parameter e(= eoi), 
or more correctly x(= Xoi)- However, e is a microscopic system parameter, independent of the 
thermodynamic state of the system; thus, it must be independent of the composition, the degree 
of polymerization, etc. On the other hand, c is determined by the thermodynamic state. Thus, 
its value will change with the thermodynamic state, the degree of polymerization, etc. Only in the 
RMA limit, see H28() below, does one find a trivial proportionality relationship between the two 
quantities c^ and e, the constant of proportionality being determined by the square of the pure 
component monomer density (f)^; otherwise, there is no extra dependence on temperature and 
pressure of the system. This RMA behavior is certainly not going to be observed in real systems, 
where we expect a complex relationship between c and e. A similar criticism also applies to the 
behavior of ci2, for which one considers the energy of mixing AEm] see below. 

1. Solubility 

In the RMA limit of an incompressible binary mixture, it is found that ci2 is proportional to the 
mutual interaction energy ei2 between the two components, see ()39() below. Thus, the sign of ci2 is 
determined by ei2. On the other hand, AEm in the RMA limit is proportional to the excess energy 

10 



e(= £12) or xi= X12)) so that its sign is determined by e. However, the solubihty of component 
1 in 2 in a given state is determined not by their mutual interaction energy ei2, which is usually 
attractive, but by the sign of the excess energy e, and the entropy of mixing. For an incompressible 
binary mixture, we have only one exchange energy e. Even away from the RMA limit for the 
incompressible mixture, a positive e implies that the two components will certainly phase separate 
at low enough temperatures. Their high solubility (at constant volume) at very high temperatures 
is mostly due to the entropy of mixing, but the energy of mixing will play an important role at 
intermediate temperatures. The solubility increases as e decreases. It also increases as T increases, 
unless one encounters a lower critical solution temperature (LCST) in which case the solubility will 
decrease with T. It is well known that LCST can occur in a blend due to compressibilty; see, for 
example, [l4|. Similarly, the solubility at constant pressure will usually decrease with temperature. 
However, it is also possible for isobaric solubility to first increase and then decrease with T. Thus, a 
properly defined mutual cohesive energy for a compressible blend should be able to capture such a 
features. A negative e implies that the two components will never phase separate. Thus, a complete 
thermodynamic consideration is needed for a comprehensive study of solubility even when we are 
not in the RMA limit, and requires investigating thermodynamic quantities such as C12 or the 
(energetic) effective chi Q to which we now turn. This is even more so important when we need 
to account for compressibility. 

2. Energetic Effective Chi 

We have recently investigated a similar issue by considering the behavior of the effective chi in 
polymers [15] ! where an effective chi, relevant in scattering experiments, was defined in terms of the 
excess second derivative of the free energy with respect to some reference state. This investigation 
is different in spirit from other investigations in the literature, such as the one carried out by Wolf 
and coworkers |l(i | where one only considers the free energy of mixing. As said above, AE'm or 
X^Q is determined by the excess energy £12, and not by eu- However, since C12 is supposed to be a 
measure of eu, it is defined indirectly by that part of the energy of mixing AEm or x^s^ ^^^ ©> 
that is supposedly determined by the mutual energy of interaction ei2. For this, one must subtract 
the contributions of the pure components from AEm ov y^g- see below for clarity. Because of this, 
even though the previous investigation of the effective chi [13] provides a clue to what might be 
expected as far as the complex behavior of x^s is involved, a separate investigation is required for 
the behavior of the cohesive density C12, which we undertake here. We will borrow some of the 
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Q, 



ideas developed in [15|, especially the requirements that 

(i) the cohesive energy density Cij for an i-j mixture vanish with Cij, and 

(ii) the formulas to determine c,-, reduce to the standard RMA form under the RMA limit. 

n 

The first condition replaces the requirement in [15|| that the effective chi vanish with £12. The 
second requirement is the same as in _15]. Its importance lies in the simple fact that any thermo- 
dynamically consistent theory must reduce to the same unique theory in the RMA limit; see |7l.ll5l| 
for details. We also borrows the idea of reference states that would be fully explained below. 

3. Symmetric and Asymmetric Blends 

It has been shown in [1^ that the non-isometric regular solution theory is more successful than 
RST but again only for a symmetric blend. A symmetric blend is one in which not only the two 
polymers have the same degree of polymerization (Mi = M2), but they also have identical inter- 
actions in their pure states (en = 622)- For asymmetric blends (blends that are not symmetric), 
even the non-isometric regular solution theory is qualitatively wrong. The si gnif icant conclusion of 
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18, and references 



the previous work is that the recursive lattice theory developed in our group 
therein], and briefly discussed in the next section to help the reader, is successful in explaining 
several experimental results where the non-isometric regular solution theory fails. A similar con- 
clusion is also arrived at when we apply our recursive theory to study the behavior of effective chi 
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1^ 



It was shown a while back that the recursive lattice theory is more reliable than the conventional 
mean field theory [l9|; the latter is formulated by exploiting RMA, which is what the regular 
solution theories are. The recursive lattice theory goes beyond RMA and is successful in explaining 
several experimental observations that could not be explained by the mean field theory. Our aim 
here is to apply the recursive theory to study solubility and to see the possible modifications due 
to 

1. finite q, 

2. non-weak interactions (e > 0), 

3. non-isometric mixing, and 

4. disparity in size (asymmetry in the degree of polymerization) and/or pure component inter- 
actions. 
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II. RECURSIVE LATTICE THEORY 

A. Lattice Model 

We will consider a lattice model for a multicomponent polymer mixture in the following, in which 
only nearest-neighbor interactions are permitted. As above, q denotes the lattice coordination 
number. The number of lattice sites will be denoted by N, and the lattice cell volume by vq, 
so that the lattice volume V = Nvq. The need for using the same coordination number and cell 
volume for the mixture and for the pure components has been already discussed elsewhere |9| in 
order to have a consistent thermodynamics. The monomers and voids are allowed to reside on 
lattice sites. The excluded volume restrictions are accounted for as described above. As shown 
elsewhere |S], one only needs to consider excess energies of interaction Sij between monomers of 
two distinct components i, and j; here, eij represents the interaction energy between two nearest- 
neighbor monomers of components i, and j to investigate the model. To model free volume, one of 
the components will represent voids or holes, always to be denoted by i = here. Thus, for a pure 
component of, for example, component i = 1, the excess interaction energy e is eoi = — (l/2)eii. 
Usually, eii is negative, which makes e positive for a pure component. Let N^i denote the number 
of monomers belonging to the ith species and N^ the number of monomers of all species, so that 
the number of voids Nq is given by A'^o = N — N^. Similarly, let Nij denote the number of nearest- 
neighbor contacts between monomers of components i, and j. The densities in the following are 
defined with respect to the number of sites. 

B. Recursive Theory 

In the present work, we will use for calculation the results developed by our group in which 
we solve the lattice model of a multicomponent polymer mixture by replacing the original lattice 
by a Bethe lattice 
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18l |. and solving it exactly using the recursive technique, which is standard 
by now. The calculation is done in the grand canonical ensemble. Thus, the volume V is taken 
to be fixed. We will assume that all material components (i > 0) are linear polymers in nature. 
The degree of polymerization of the zth component, i.e. the number of consecutive sites occupied 
by it, is denoted by Mj > 1. The linear polymers also include monomers for which Mi = 1. Each 
void {i = 0) occupies a single site on the lattice. Let 0o = Nq/N denote the density of voids, 
<Pmi = N^i/N the monomer density of the zth component, and (pij = Nij/N the density of nearest- 
neighbor (chemically unbonded) contacts between the monomers of the two components i, and j. 

13 



It is obvious that 

4>ia = '^ (prai = 'i- - (1)0 (16) 

j>0 

denotes the density of ah material components (i > 0) monomers. The density of all chemical 
bonds is given by 

(l> = ^(l)mm, yi = {Mi-l)/Mi. (17) 

i>0 

The quantity 0u = (7/2 — (j) denotes the density of lattice bonds not covered by the polymers. Let 
us also introduce Qi = q — 2vi, (f>iu = gi(/)mi/2. As shown in [a, |l^, the pressure P is given by 

P = {kBT/vQ)[-\n4>Q + {q/2)\u{2<i)Jq)] + (A:Brg/2t>o) ln(</.[Jo/</.oo)], (18) 

where k-Q is the Boltzmann constant, (J)qq is the density of nearest-neighbor void-void contacts, and 
(^QQ is its athermal value when all excess interactions e^j are identically zero. The athermal values 
of (j)ij are given by 

4>^j = 2(f)iu(t)ju/4'u{'^ + Sij), 

where 5ij is the Kronecker delta, and give the values of the contact densities in the atermal state 
when all Eij = 0. The athermal state has the same volume V as of the original system. The first 
term in ()18|) give the athermal value Path of P, and the second term is the correction 

Pint = {kBTq/2vo) HcP^oo/cPoo) (19) 

to the athermal pressure due to interactions. For attractive interactions responsible for cohesion, 
Pint is going to be negative, as discussed above. This is the correction to Path due to interactions and 
determines the internal pressure Pin = —Pint- The identification also holds for pure components, 
so that Pin remains positive for attractive interactions. 

Since there is no kinetic energy on a lattice, the internal energy in the lattice model is purely 
due to interactions, and this energy per unit volume -Eint is given by 

-E^int = ^ eij(j)ij/vo, (20) 

i>j>0 

which will be used here to calculate the cohesive energy density, also known as the cohesive pressure. 
Note that the form of the pressure in (jlHj) does not explicitly depend on the number of components 
in the mixture. It should also be clear from (|18|) that the incompressible state (P -^ oo) corresponds 
to (/>o ^ at any finite temperature. However, we will not be interested in this limit in this work. 
In all our calculations, we will keep (/>o > 0. 
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C. RMA Limit 

The approximation has been discussed in details in |9|,|l3|, so we wih only summarize the results. 
This limit is very important since all thermodynamically consistent lattice theories must reduce 
to the same unique theory in the RMA limit. Thus, the RMA limit provides a unique theory, and 
can serve as a testing ground for the consistency of any theory. We now show that our recursive 
theory reproduces the known results in the RMA limit. To derive this unique theory, we note that 
in this limit, we have qi > q, and that the contact densities take the limiting form 

'P'ij * q(prai4'mj/il+Sij),4'ij ' q4'iai4>mj / {I + Sij) , q ln{2(j)^ / q) * - 2(p. (21) 

Finally, using (pj^^Q to also denote the void density (pQ, we have in this limit 

Eint ' ^ qeij(prai4>iaj/'^Vo, (22a) 

PPvo • -Incpo- (j)-'^X0i4'rai+ ^ Xij4'rai(pmj- (22b) 

i>0 j>i>0 

For the pure component (ith component) quantities, we use an additional superscript (P) for 
some quantities, such as (pf^ representing the contact density between unbonded monomers, or use 
the superscript i for_pther quantities, such as E^J^ representing the pure component internal energy. 



In the RMA hmit \9\, we obtain 

njjiii) RMA iP2 / ar,{i) RMA ^P2 / /oo\ 

by restricting ()22al22b|) to a single material component i > in addition to the species 0. Here, we 
have used the fact that xoi = — (l/2)/3'ZeM, and that 1 — cp^^ represents the pure component free 
volume density. We notice that in the RMA limit, 

p(i) RMA p(i) 
-^int Mnt 

for a pure component. But this equality will not hold when we go beyond RMA. 

D. Infinite Temperature Behavior 

In the limit /3 — > at fixed e^j , the limiting form of E'int , and P are 

i>j 

l3Pvo ^ pP^thVo = -\ncl>o + {q/2) ln(2</.,/g). 
15 



and 

For any finite pressure P , we immediately note tliat (/)o — > 1, that is the entire lattice is covered by 
voids with probability 1. This shows that at a fixed and finite pressure P, (p and 0^-, and therefore, 
E'int vanish as T ^ oo. On the other hand, if the volume is kept fixed, which requires the free 
volume density (po to be strictly less than 1, then P -^ oo, as T ^ oo, and cj) and </>^-, and therefore, 
Sint does not vanish as T ^ oo. Thus, whether there is cohesion or not at infinite temperatures 
depend on the process carried out. This will emerge in our calculation below. 

E. Choice of parameters for Numerical Results 

In the following, we will take vq = 1.6 x 10^^^ m^ for all calculations. We will take en = 
—2.6 X 10~^^ J, and 622 = —2.2 x 10^^^ J for the two components, unless specified otherwise. The 
degree of polymerization M will be allowed to take various values between 10, and 500, and q is 
allowed to take various values between 6, and 14 as specified case by case In some of the results, we 
will keep the product eq fixed, for example eq = 1.56 x lO"^'^ J, as we change q, so that e changes 
with q, since the cohesive energy density is determined by this product; see below. 

III. INTERNAL PRESSURE AND (d£/dV)T 

The internal pressure Pin = —-Pint is obtained by subtracting Pathi see (fTl)) . from the isentropic 
volume derivative of the energy P = —(d£/dV)s, which follows from the first law of thermody- 
namics; the derivative is also supposed be carried out at fixed number of particles A'p. We will 
assume here that Path is the pressure of the athermal system at the same volume as the real system 
(Kth = V). 

The internal pressure is not the same as the isothermal derivative —{d£/dV)T- To demonstrate 
this, we start with the thermodynamic identity 

-{d£/dV)T = P - T{dP/dT)v. 

The derivatives here and below are also at fix A^p. We now note that the internal energy is a 
sum of kinetic and potential (interaction) energies K and [/, respectively. If U does not depend 
on particle momenta, which is the case we consider here, then the canonical partition function 
in classical statistical mechanics becomes a product of two terms, one of which depends only on 
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particle momenta, and the other one on U. Thus, the entropy S of the system in classical statistical 
mechanics is a sum of two terms 5ke and Sconi, where Ske depends on K, and 5conf on U; see for 
example, J20|. The entropy 5ke due to the kinetic energy depends on K and is independent of the 
volume and the interaction energy. The configurational entropy 5conf of the system, on the other 
hand, is determined by U and is a function only of the volume V when there is no interaction 
([/ = 0), i.e. in the athermal state. Thus, for a system in the athermal state, 

S^K, V) = 5eo„f,ath(^) + Ske{K), 

where K{T) is the kinetic energy of the system at a given temperature T [20|, so that 

{dS^th/dV)K = {dSconiMh/dV) = /?Path. (24) 

Since the identity (|24|) is valid in general for the athermal entropy, we conclude that /3-Path is a 
function independent of the temperature. Rather, it is a function of A'^p, and V; indeed, it must 
be simply a function of the number density np = Np/V of particles. Consequently, the athermal 
pressure Path depends linearly on the temperature, so that 

T{dP^th/dT)v = Path. 
Using this observation, we find that 

{d£/dV)T = Pin - T{dPm/dT)v = [^(/3PiN)]y, (25) 

clearly establishing that Pin / {d£/dV)T, unless {dPi^/dT)v vanishes, which will happen if Pint 
becomes independent of the temperature (at constant volume). (As we will see below, this will 
happen in the RMA limit.) It is also possible for Pn and {d£/dV)T to be the same at isolated 
points, the extrema of Pin as a function of T, where {dPi-^/dT)y = 0. 
In the RMA limit, we have 

Pin = -Pmt ^^ - qe<f'j2vo, (26) 

so that {dP\at/dT)v = 0. Thus, we conclude that 

Pin '^ {d£/dV)T, and c^ ^^ Pin 

in the RMA limit, as claimed earlier. 

Considering the RMA behavior of Pm, and its equality with c^ in the same limit, we infer that 
we can also use the internal pressure Pin as an alternative quantity to measure cohesiveness, as 
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was first noted by Hildebrand [2]- However, in general, the two quantities Pin and c are not going 
to be the same. 

We can directly calculate the internal pressure and its temperature-dependence either at con- 
stant volume or at constant pressure in our theory from H18|). The internal pressure Pin is given by 
the negative of Pint in (|19jl . and can be used as an alternative quantity that is just as good a mea- 
sure of the cohesion as the cohesive energy density [2], as was discussed above. In Fig. ^ we show 
the interaction pressure Pint^o for a symmetric blend (M = 100) as a function of kT. Both axes 
are in arbitrary energy unit. In the same energy unit, the excess energies are en = —0.25 = 622, 
and ei2 = —0.249, and the coordination number is g = 6. The free volume density is kept fixed 
at 00 = 0.1, as we change the temperature. This means that we also keep the total monomer 
density fixed. Thus, if we consider the case of a fixed amount of both polymers, then this analysis 
corresponds to keeping the volume fixed as the temperature is varied. In other words. Pint in Fig. 
n is for an isochoric process. We immediately notice that Pint not only changes with T, so that 
the second term in ()25() does not vanish, but it is also non-monotonic. Moreover, the difference 
between {d£ /dV)T and Pin could be substantial, especially at low temperatures (see Fig. P) so 
that using {d£/dV)T for Pin could be quite misleading. 

The minimum in Pint occurs at kT ~ 2.0 in Fig. ^ At this point, {d£/dV)T and Pin become 
identical, but nowhere else. However, they become asymptotically close to each other at very high 
temperatures where the slope {dPm / dT)v gradually vanishes. Let us fix the arbitrary energy 
unit to be 1.38 x lO^^^J (equal in magnitude to lOO/c), so that kT = 2.0 in this energy unit, the 
approximate location of the minimum in Pint in Fig. ^ The minimum occurs around T = 200 
K, i.e., tQ = —73.15 °C. Thus, Pin is an decreasing function of temperature approximately above 
tc — —73.15 °C, and an increasing function below it. The location of the maximum in Pin will 
change with the applied pressure, the two degrees of polymerization, the interaction energies, etc. 
The significant feature of Fig. ^ is the very broad flat region near the minimum, which implies a 
very broad flat range in temperature tc near the maximum of Pin . 

In Fig. 12 we show Pin (MPa) as a function of the temperature tc (°C) for a pure component 
(M = 100, q = 10, vq = 1.6 X lO'^s m^, and e = -2.6 x lO^^i j) foj. ^ constant V (0o = 0.012246), 
and constant P (1.0 atm). At 0°C, the system has 0o = 0.012246, and P = 1.0 atm in both 
cases so that isobaric and isochoric Pin match there, as seen in Fig. \^ We notice that isochoric 
Pn has a very week dependence on tc (suggesting that the temperature range in Fig. [2lmay be 
near the maximum of Pin)) while isobaric calculation provides a strongly dependent Pin, which 
asymptotically goes to zero. This difference in the isobaric-isochoric behavior is consistent with 
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FIG. 1: Non-monotonic behavior of uoPint with kT (arbitrary energy scale). 

our claim above. The difference between the two internal pressures finally reaches a constant equal 
to about 80 MPa at very high temperatures. We also show the derivative {d£/dV)T calculated for 
the isochoric case for comparison. We use isochoric Pin and (|25j) to obtain isochoric {d£/dV)T- 
We notice that it differs from isochoric Pin by a small amount over the entire temperature range 
considered. Near 0°C, they differ by about 1 MPa, and this difference decreases as the temperature 
rises so that they approach each other at higher temperatures. This is consistent with what we 
learned from Fig. ^ above. We also notice that Pin > {d£/dV)T over the temperature range in 
Fig. [21 and the difference gradually vanishes. From (^S)). we conclude that {dPm/dT) > for this 
to be true. This corresponds to (dPint/dT) < in Fig. |21 Thus, the temperature range is below 
the minimum of Pint! refer to Fig. ^ 

IV. PURE COMPONENT COHESIVE ENERGY DENSITY OR PRESSURE 
A. Athermal Reference State 

The following discussion in this section is general and not restricted to a lattice except for 
the numerical results, which are on a lattice. Therefore, we will include the kinetic energy of the 
system in this discussion. The definition of 5 requires considering a pure component consisting of 
one particular component. To be specific, we consider this to be the component i = 1, and we will 
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FIG. 2: Various isochoric and isobaric pressures vs temperature for a pure component (Af — 100, q — 10, vq = 
1.6 X 10^^*m'^,and e = —2.6 x 10^^^ J). In the inset, we show cohesive pressure at 1.0 atm for two different 
M (Af = 100 :componcnt 1 and M = 10 :componcnt 2), with the smaller M showing a discontinuity at its 
boiling point around 600° C. 



not exhibit the index unless needed for clarity. At a given temperature T, and pressure P or volume 
V, one calculates the interaction energy ii^int per unit volume by subtracting the kinetic energy K 
due to motion from the total internal energy £ of the system. We will assume in the following that 
the interaction energy depends only on the coordinates and not on the momenta of the particles, 
n classical statistical mechanics, K is a function only of the temperature T but not of the volume 
20|, and can be obtained by considering a fictitious pure component at the same T, and P (or 
V), containing the same number of particles, but with particles having no interaction (U = 0). 
This fictitious reference system is, as said above, known as the athermal reference state. Then, 
K = <Sath) where i?ath is the energy of the athermal system. This reference state is equivalent to the 
conventional view according to which one considers the gaseous state of the system at the same 
T, but at zero pressure, so that y ^ cxd, which allows for infinite separation between particles to 
ensure Sint = 0. The energy of the gaseous state is exactly <?ath; which is independent of the volume, 
even though the physical state requires ^ — > oo for the absence of interactions. We will continue 
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to use the athermal system instead of the gas phase since the latter requires considering different 
pressure or volume than the pressure or volume of the physical system under consideration. Thus, 
we find from (^ 

Eint = E — iys_th/y)Ea.thi 

where E = 8 /V is the energy density of the pure component (with interactions), and -Eath = 
"^ath/Kith is the energy density of the athermal reference state (e = en = 0), and V and T4th 
are the corresponding volumes. As said above, -Emt has different limits at infinite temperatures 
depending on whether we consider an isochoric or isobaric process. 

B. Cohesive Energy Density 

The cohesive density c^ is the interaction energy of the particles in the pure component. We 
set e = en, which must be negative for cohesion. The cohesive energy density c for the pure 
component can be thought of as functions of T, and P, or T, and V, as the case may be. It is 
obvious that i^int vanishes as the interaction energy e vanishes. This will also be a requirement for 
the cohesive energy density: c should vanish with the interaction energy. We have calculated c 
as a function of T for isochoric (constant volume) and for isobaric (constant pressure) processes. 
We denote the two quantities by (^ and (^, respectively, and show them in Fig. |21 where they 
can also be compared with Pjn for isochoric and isobaric processes, respectively. The conditions 
for the two processes are such that they correspond to identical states at 0°C. We find that while 
the two quantities c and Pin are very different for the two processes, they are similar for each 
process. We again note that Cy and Cp behave very differently, as discussed above, with (^ > (^. 
Not surprisingly, the same inequality also holds for Pn- 

The almost constancy of isochoric c^ and Pin provides a strong argument in support of their 
usefulness as a suitable candidate for the cohesive pressure. This should not be taken to imply 
that (? and Pin remain almost constant in every process, as the isobaric results above clearly 
establish. Unfortunately, most of the experiments are done under isobaric conditions; hence the 
use of isochoric cohesive quantities may not be useful, and even misleading and care has to be 
exercised. 

In the lattice model with only nearest-neighbor interactions, -Eint = &<Pc where (j)c = <pii is the 
contact density between monomers (of component i = 1); hence 

6'^ = c^ = -e(l)JvQ. (27) 
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a. RMA Limit In the RMA limit, we find from ((231) that 

c^m^ - ge</2.o, (28) 

where (jj-ca represents the pure component monomer density in this section. We thus find that the 
ratio 

c Iqecht, = const 



in this hmit. Since in this hmit, the product qe remains a constant, the ratio 

c = c /cf)^ = const 

in this hmit. Both these properties wih not hold true in general. 

h. End Group Effects In the following, we will study the combinations 

c= 2c^vo/iq - 2i^)e</'m,c = c^ /{q - 2v) (29) 

which incorporate the end group effects along with the linear connectivity via the correction 2z/, 
where 

v = {M - 1)/M. 

In the RMA hmit, 

C > 1, 

and 

/ , RMA ,2 /o 

c/e ► </'m/2vo- 



We will investigate how close c is to 1 in our calculation for finite q, and how close to a quadratic 
form in terms of (jy^a does c' have. 

C. van der Waals Fluid 

To focus our attention, let us consider a fiuid which is described by the van der Waals equation. 
It is known that for this fiuid, -Emt = —Npu/V'^, where a > is determined by the integral |2ll | 

a = -(1/2) / udV; (30) 
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here u is a two-body potential function. It is clear that a is independent of T for the van der Waals 
fluid. We can also treat it to be independent of V (any ^-dependence must be very weak and can 
be neglected). The lower limit of the integral is the zero of the two-body potential energy. Thus, 

c = arip, (31) 

where np is the number density per unit volume; for polymers, np = (f)^/MvQ. It is interesting to 
compare (|31j) with (|28|) derived in the RMA limit. They both show the same quadratic dependence 
on the number density. We also note that c^ vanishes as a vanishes, as we expect, but most 
importantly the ratio c^/np is independent of the temperature, just as c is a constant in the RMA 
limit. Deviation from the quadratic dependence will be observed in most realistic systems, since 
they cannot be described by this approximation. 

An interesting observation about this fluid is that Pint is exactly equal to -Emt! thus, (? = —Pint 
for this fluid. In addition, it is well known that for this fluid 

{d£/dV)T = {d£,^t/dV)T = anl, (32) 

since a is temperature-independent. Thus, 

Pin = {d£/dV)T (33) 

for the the van der Waals fluid, in which a must be taken as T- independent. The equality ()33|) is 
not always valid as discussed above. 

D. The Usual Approximation 

Usually, one approximates E'int by the energy density of vaporization at the boiling point at 
r = Tb . This means that one approximates c^ (T, P) by its value c^ (Tb , P) , which will be a function 
of P, but not of T. On the contrary, c {T,P) will show variation with respect to both variables. 
In addition, it will also change with the lattice coordination number q, and the interaction energy 
e. It is clear that isobaric c (T,P) will show the discontinuity at the boiling point, as shown in 
the inset in Fig. [21 where we report isobaric c for two different molecular weights M = 10, and 
M = 100 at 1.0 atm. The shorter polymer system boils at about 600°C (and 1.0 atm), but not 
the longer polymer over the temperature range shown there. At 0°C, the smaller polymer has a 
slightly higher value of isobaric c^. If we calculate isochoric c^ at a volume equal to that in the 
isobaric case at 0°C and 1.0 atm, then the corresponding isochoric c , which is almost a constant 
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FIG. 3: Cohesive energy density as a function of pressure for different q. We also show c' and c. 

(as shown in the main figure in Fig. EJ, will be about 65 MPa, close to its value at 0°C. This 
value is much larger than c (Tb, P) of about 15 MPa at the boiling point. Thus, the conventional 
approximation cannot be taken as a very good estimate of the cohesiveness of the system, which 
clearly depends on the state of the system. Cohesive energy density as a function of pressure for 
different q. 

E. Numerical results 

We will evaluate various quantities for isochoric and isobaric processes. For isochoric calcula- 
tions, we calculate 0m at some reference T,P. Usually, we take it to be at 1°C, and 1.0 atm. We 
then keep (j)^ fixed as we change the temperature. This is equivalent to keeping the volume fixed 
for a given amount of the polymers. For isobaric (constant P) calculations, we again start at 1°C, 
and take a certain pressure such as 1.0 atm, and keep the pressure fixed by adjusting 0m as we 
change the temperature. For a given amount of polymers, this amounts to adjusting the volume 
of the system. 

We show c^ in Fig. |31 and the monomer density 0m in Fig- El as a function of P (P < 2 atm) 
for different values of q from g = 6 to g = 14. We have taken M = 500, vq = 1.6 x 10~^^ m^, and 
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FIG. 4: Monomer density as a function of presssure for different q. For q > 8, we have a liquid which gets 
denser as q increases. For q — Q^ we have a gas phase at low pressure which becomes a lighter liquid at 
higher pressure, with a liquid gas transition around 0.05 atm. 

e = —2.6 X 10^^^ J, and have set tc = 500°C. We see that (? increases with q, which is expected; 
see ()28() . To further analyze this dependence, we plot c' in Fig. 3. There is still a residual increase 
with q. The increase is also partially due to the fact that (j)^ increases with q, as shown in Fig. 
|1I Therefore, we also plot c as a function of P in Fig. |21 We notice that there is still a residual 
dependence on q with c increasing with g, and reaching 1.0 from below as q increases, which is 
consistent the RMA limit. In Fig. 13 we plot c and (/)m (in the inset) for the same system over a 
much wider range of pressure for different q. From the behavior of c and ^m, we easily conclude 
that (? changes strongly with P for smaller q, and the dependence gets weaker as q increases. 

In Fig. IHl we plot (? for q = 12 as, a. function of the monomer density (f)^. We have taken 
M = 500, vq = 1.6 X 10"28 m^, and e = -2.6 x lO^^i j^ ^nd have set tc = 500°C. According to 
H28|) , it should be a quadratic function of the monomer density. To see if this is true for the present 
case of finite q, we plot the ratio c /(j)'^ in the inset, which clearly shows that there is still a strong 
residual dependence left in the ratio. Thus, c^ in a realistic system should not be quadratic in (f)^. 
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FIG. 5: c as a function of pressure at 500°C for different q. As q increases, c approaches 1, its limiting RMA 
value, which also forms its upper limit. For any finite g, c is strictly lower. 

In Fig. [3 we show c as a function P for small pressures for different M ranging from M = 10 
to M = 500; we keep q = 14, vq = 1.6 x 10^^^ m'^, and e = —2.6 x 10^^^ J, and have set the 
temperature tc = 25° C. We see that c^ decreases as M increases, but different curves are almost 
parallel, indicating that it is not the slope, but the magnitude that is affected by M. In Fig. |H1 
we show c as a function P for small pressures for different temperatures ranging from tc = 25° C 
to tc = 5000°C; we keep g = 14, M = 100, vq = 1.6 x lO^^s ^^3^ ^^^ ^ ^ _2q x lO^^i j_ ^g 
immediately note that the pressure variation is quite minimal over such a small range of pressure 
between 1 to 25 atm. However, the temperature variation is quite pronounced, again affecting 
the magnitude and not the slope. There are two important observations, (i) For the highest 
temperature 5000°C, the state corresponds to a gas, since (? ~ 0. (ii) At low temperatures, c^ 
reaches an asymptotic value since the system has become almost incompressible. 

In Fig. 13 which is for a pure component studied in Fig. |21 we keep the product e{q/2 — u) fixed 
at (—5.21 X 10~^-'^J), as we change q, so that e changes with q. We keep M = 500, vq = 1.6 x 10"^^ 
m^, and tc = 500°C. The effect of changing q is now minimal. The cohesive energy density shows 
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FIG. 6: c^ as a function of (f)j-n. In the inset, we show c, showing that c^ is not quadratic in (/),n, as suggested 
by the regular solution theory. 

minimal variation with q, except that it is higher for larger q, but the difference first increases with 
(pm and then decreases as 0m — > 1. In Fig. I10( we keep instead the product eq fixed, so that there 
is no endpoint correction. In contrast with Fig. |^ we find that c^, while still increasing with q, 
has the property that the its difference for different q continues to increases with <p^ as ^m -^ 1. 
It is clear from both figures that c^ increases with g at a given (/)m and saturates as q becomes 
large, which is the direction in which q must increase to obtain the RMA limit. Thus, c achieves 
its maximum value at a given density </)ni in the RMA approximation: for any realistic system in 
which q is some finite value, the value of the cohesive energy at that density </>m will be strictly 
smaller. It is also evident from the inset in both figures that at a given pressure P, c^ achieves its 
minimum value in the RMA approximation. 

The results for isochoric and isobaric calculations are shown in Fig. 1111 We consider q = 14, 
M = 100, vq = 1.6 X 10^^*^ m'^, and e = —2.6 x 10~^^ J, and have taken the pressure to be P = 1.0 
atm at the initial temperature tc = 0°C. This corresponds to the monomer density 0m = 0.997115. 
For the isochoric calculation, we maintain the same monomer density (0m = 0.997115) at all 
temperatures. For the isobaric calculation, we maintain the same pressure (P = 1 atm) at all 
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FIG. 7: Cohesive energy density as a function of pressure for different molecular weights. 
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FIG. 8: Cohesive energy density as a function of pressure at different temperatures. 
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FIG. 9: c^ as a function of (p^i for different q when e{q/2 — i') is kept fixed at (—5.21 x 10 ^^J). 

temperatures. This ensures that isochoric Cy and isobaric c^ are identical at tc = 0°C. The figure 
clearly shows that Cy is always higher than Cp, as discussed above. This behavior is not hard to 
understand. As we raise T, the pressure increases in the constant volume calculation above P = 1 
atm, making particles get closer together, thereby increasing the cohesive energy. 

V. MUTUAL COHESIVE ENERGY DENSITY OR PRESSURE 
A. Mutual Interaction 

The solubility of one of the components in a binary mixture of components i = 1,2 increases 
as the excess energy ei2 decreases for the simple reason that a larger ei2 corresponds to a stronger 
excess repulsion between the two components. In particular, the two components will not phase 
separate at any temperature if £12 < 0. If £12 > 0, the two components will phase separate at low 
temperatures. Thus, knowing whether £12 > immediately allows us to conclude that solubility 
will not occur everywhere. 

Usually, all energies eij are negative, but the sign of £12 depends on their relative magnitudes. 
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FIG. 10: c^ as a function of 0,n for different q when eq is kept fixed at (—1.56 x 10 ^°J). 



If, however, we assume the London conjecture (|11|) . then the corresponding excess energy becomes 



£12 = iV\e^\ - v1ed)V2 > 0. 



(34) 



Thus, the London conjecture implies that the two components experience a repulsive excess inter- 
action so that the solubility decreases as £12 increases, i.e., as en and 622 become more disparate. 
The maximum solubility in this case occurs when £12 = 0, i.e., when the two components have 
identical interactions ("like dissolves like"). In this case, the size or architectural disparity cannot 
diminish the solubility because the entropy of mixing will always promote miscibility. In general, 
£12 > 0, and it becomes necessary to study solubility under different thermodynamic conditions. 

To simplify our investigation, we will assume the London conjecture (|llj) for the mixture in all 
the calculations. Thus, the two components always experience a repulsive excess interaction in our 
calculation. (This is also true when we intentionally violate the London conjecture ()11() and set 
ei2 = 0, as is the case for SRS.) In the RMA limit, the conjecture (fTT|) immediately leads to (fTn|l. 
as we have seen above. This need not remain true when we go beyond RMA. Thus, we will inquire 
if ()lUp is satisfied in general for cohesive energies that are calculated from our theory under the 
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FIG. 11: Various c^ as a function of temperature. 

assumption that the London conjecture (|11() is vahd. Any failure of (jlUp under this condition wih 
clearly have significant implications for our basic understanding of solubility. 

B. van Laar-Hildebrand Approach using Energy of Mixing 

We follow van Laar l] and Hildebrand 3], and introduce Ci2 by exploiting the energy of mixing 
A£^M per unit volume. According to the isometric regular solution theory [3], the two are related 
by 



A£^M = (cn + C22 - 2ci2)v?i¥'2, 



(35) 



where ipi are supposed to denote the volume fractions. Using the London-Berthelot conjecture (|1()|) . 
we immediately retrieve (jSJ, once we recognize that df = Cii, and 82 = C22, see ^ above. We now 
take ()35() as the general definition of the mutual cohesive energy density C12 in terms of the pure 
component cohesive energy densities. The extension then allows us to evaluate C12 by calculating 
AEm, provided we know c^i and C22 for the pure components; the latter are independent of the 
composition. It can be argued that since (|35() is valid only for isometric mixing, it should not be 
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FIG. 12: Mutual cohesive energy density as a function of temperature for different q. We have taken a blend 
at 50-50 composition (Mi = Ah = 100; en = -2.2 x lO^^i j^ ^^^ ^ _2q x IO-21 j^ y^ = i.g x IQ-^s ^3) 
at 1.0 atm. 

considered a general definition of C12 for non-isometric mixing. However, since one of our objectives 
is to investigate the effects due to isometric and non-isometric mixing, we will adopt H35|) as the 
general definition of ci2- 

1. Lattice model 

The kinetic energy of the mixture is the sum of the kinetic energies of the pure components, all 
having the same temperature, and will not affect the energy of mixing. Thus, we need not consider 
the kinetic energy in our consideration anymore. In other words, we can safely use a lattice model 
in calculating AE'm- This is what we intend to do below. 

The definition of C12 depends on a form (|35|) whose validity in general is questionable, as it is 
based on RST. To appreciate this more clearly, let us find out the conditions under which the RMA 
limit of our recursive theory will reproduce ()35() . We first note that the interaction energy density 
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(per unit volume) E'int of the mixture from (|2U() is given by 

-EintfO = 611011 + e22<A22 + ei2<?!'l2 = -Cn^'O'All/'/'n - C^2''^o4'22/4>22 + 612012, (36) 

where we introduced the pure component cohesive energy densities in the last equation. The energy 
of mixing per unit volume is 

AEmVq = 611011 + 622022 + 612012 " (^l^/^)6ll0n " (V2^/V')e22022, (37) 

where V^^^, and V^^ are the pure component volumes. It is easy to see that in general 

V[/V = X0^/0Pi ^ 0ml/0!;".l, V^/V = (1 - X)0„,/0P2 ^ 0m2/0L, (38) 

where x = 0mi/0m is the monomer fraction of species i = 1 introduced earlier in @,and 0j^j the 
pure component monomer density of the ith species. The monomer density of both species in the 
mixture is 0m = 0ml + 0m2- 

2. RMA Limit and Monomer Density Equality: 

In the RMA limit, it is easy to see by the use of (|22b|) that the last equation in (|36|) reduces to 

^.nt ^ -cUvr/V?-cl,{Vi'/Vf-2cMVr/V){Vr/V), 

where 

cfi > - {l/2)qeii4>^i/vo, cu ' - qeu4>mi4>m2/'^vo (39) 

are the cohesive energy densities in this limit. The form of cj^ is exactly the same as the one derived 
above for the pure component; see H23|). and H31|). It is also a trivial exercise to see that the RMA 
limit of the energy of mixing 1)3 7|) will exactly reproduce ()35() provided we identify 

and further assume 

vi + vi = V, 

so that <pi+ ip2 = 1- The last condition is nothing but the requirement that mixing be isometric 
and can be rewritten using (|38j) as 

^m0ml + a;0m(0m2 " (Ami) = 0ml (^L- 
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This should be vahd for all x including x = and x = 1. This can only be true if we require the 
monomer density equality: 

(t>ra = '/'ml = </'m2- (40) 

This condition is nothing but the equality of the free volume densities in the mixture and the pure 
components, and is a consequence of isometric mixing, as is easily seen from (|45a() obtained below. 
We finally conclude that 

ifl=X,if2 = l-X, (41) 

as was also the case discussed earlier in the context of ©• 

3. Isometric RMA Limit 

Let us now consider the isometric RMA limit for which we have the simplification 

Vl/V = X, V^/V = l-x- 
see (001. Thus, 

j;, RMA r P 2 , o I P 2i 

as is well known 3]. This form can only be justified in the isometric RMA limit, with the volume 
fractions given by (|4H) . Similarly, 

AEm-^O ' q[ei2 - (l/2)(eii + e22)]<?iml</'m2 = (x//3)^ml<?^m2, 

where we have introduced the Flory-Huggins chi parameter x = qP^- Using (|4U() . we can rewrite 
the above energy of mixing in the form ()35() . We finally have for Aii^M in the isometric RMA limit 

again with the volume fractions given by (|41() . 

4-. Volume Fractions 

In the following, we will always take ipi to be given by H41|). This ensures that 9J1+ (^2 = 1- 
Another possibility is to define (pi in terms of partial monomer volumes vl: 

Vi = (pmiVi/vo. (43) 
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However, as shown in [1^, the error is not significant except near x = or x = 1. Since the 
calculation of UJ is somewhat tedious, we will continue to use (|4ip for ipi in our calculation, as we 
are mostly interested in x = 0.5. 

5. Beyond Isometric RMA 

Beyond isometric RMA, the energy of mixing will not have the above form in (|42j) : rather, it 
will be related to the energetic effective chi introduced in @ [a, Isl in exactly the same form as 
above: 

AEm = {x^s(l)l,/Pvo)ipiip2, 

which ties the concept of cohesive energy density intimately with that of the effective chi, as noted 
earlier. However, it is also clear that cu and x^s ^^^ ^°^ directly proportional to each other. 

C. van Laar-Hildebrand ci2 

Using (P7|) for pure component cohesive energies, we obtain 

C12V0 = (612/2)012 - en [((All - V^ /V<P^^)/2^iV2 - <An] - £22 [(fe - ^2^/^</'22)/2^1^2 - 4>22] (44) 

as the general expression for the cohesive energy density. We will assume that (pi are as given in 
(|4H) . It is clear that the definition of C12 given in l|35|) is such that it not only depends on the 
state of the mixture, but also depends on pure component states. This is an unwanted feature. In 
particular, C12 will show a discontinuity if a pure component undergoes a phase change (see Fig. 
1221 later), even though the mixture does not. 

In Fig. 1121 we show C12 for a 50-50 blend (M = 100) at 1.0 atm as a function of tc calculated 
for different values of q from q = 6 to q = 14. The curvature of C12 gradually changes at low 
temperatures from concave upwards to downwards. As not only the magnitude but also the shape 
of C12 changes with q, C12 is not simply proportional to 5; it is a complicated function of it. This is 
easily seen from the values of 012/(7 ^^ ^c = 0°C, which is found to increase with q. It is about 4.0 
at q=6 and increases to about 6.4 at q = 14. The temperature where C12 asymptotically becomes 
very small, such as 0.1 MPa occurs at higher and higher values of the temperature as q increases. 
This is not surprising. We expect the cohesion to increase with q at a given temperature. 
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D. Isometric Mixing: EDIM and DDIM 

The volume of mixing is defined as 

AVM = V-V[-Vi'. 

Using (|55)) . we find that the volume of mixing per unit volume {Avu = AVu/V), and per monomer 

(AuM = AVm/Njo) are 

AVM = (pra[l/(pra " a^/<Aml " (1 " 2;)/</'m2]; (45a) 

A9m = (l/0m - x/4>^^ - (1 - x)/(l)^2)vo- (45b) 

One of the conditions for RST to be valid is that this quantity be zero (isometric mixing). The 
condition (|4Uj) . which ensures isometric mixing at all x, is much stronger than the isometric mixing 
requirement at a given fixed value of x. There are situations in which mixing is isometric at a 
given X, but ()4U() is not satisfied. For a given pure component monomer densities 4>^i, and 0^2! 
that are to be mixed at a given composition x, we must choose the mixture density (p^ to satisfy 

l/cPm = x/4>l^ + (1 - x)/4>l^ (46) 

in order to make the mixing isometric; see ()45ap . In this case, ()35() will not hold true despite mixing 
being isometric. In our calculations, we will consider both ways of ensuring isometric mixing. We 
will call the mixing method satisfying (|4U|) equal density isometric mixing (EDIM), and the mixing 
method satisfying H46|) different density isometric mixing (DDIM). 

For most of our computation, we fix the composition x. Almost all of our results are for a 50-50 
mixture. We will consider both isometric mixing processes noted above in our calculations. A 
variety of processes can be considered for each mixing. In order to make calculations feasible, we 
need to restrict the processes to a few selected ones. We have decided to investigate the following 
processes with the hope that they are sufficient to illuminate the complex behavior of cohesive 
energy densities and their usefulness. 

1. (Isometric) Isochoric Process 

The process should be properly called an isometric isochoric process, but we will use the term 
isochoric process in short in this work. The volume of the mixture is kept fixed as the temperature 
is varied. The energy of mixing can be calculated for a variety of mixing processes. We have 



36 



decided to restrict this to isometric mixing. We calculate the mixture's monomer density (jj-ca at 
each temperature. For each temperature, we use (|40p to determine the pure component monomer 
densities to ensure isometric mixing for the selected x. 

2. (Isobaric) EDIM Process 

The process should be properly called an isobaric EDIM process, but we will use the term 
EDIM process in short in this work. We keep the mixture at a fixed pressure, which is usually 1.0 
atm, and calculate its monomer density (p^^ at each temperature. We then use EDIM to ensure 
isometric mixing at each temperature and calculate the energy of mixing. In the process, the 
mixture's volume keeps changing, and the pure component pressures need not be at the mixture's 
fixed pressure. Thus, the mixing is not at constant pressure, even though the mixture's pressure 
is constant. 

3. DDIM Process 

For DDIM, we keep the pressures of the pure components fixed, which is usually 1.0 atm, and 
calculate (p^i, and 0^2' from which we calculate (p^i using (|1^ for the selected xv. This ensures 
isometric mixing, but again the mixing process is not a constant pressure one since the mixture 
need not have the same fixed pressure of the pure components. Even though the energy of mixing is 
calculated for an isometric mixing, it is neither calculated for an isobaric nor an isochoric process. 

All the above three mixing processes correspond to isometric mixing at each temperature. Thus, 
we can compare the calculated van Laar-Hildebrand ci2 for these three processes to assess the 
importance of isometric mixing on ci2- It is not easy to consider EDIM for an isochoric (constant 
V) process. Therefore, we have not investigated it in this work. 

4- Isobaric Process 

In this process, the pressure of the mixture is kept constant as the temperature is varied. The 
calculation of the energy of mixing can be carried out for a variety of mixing. We will restrict 
ourselves to mixing at constant pressure so that the pure components also have the same pressure 
as the mixture at all temperatures. Thus, the volume of mixing will not be zero in this case. The 
process is properly described as a constant pressure mixing isobaric process, but we will use the 
term isobaric process in short in this work. 
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FIG. 13: Energy of mixing and mutual cohesive energy density as a function of composition x. We take 
Ml = 10,1/2 = 100, en = 622 = ei2 = -2.6 x lO^^^ J, vq = 1.6 x lO^^s ^-^3 ^^^ ^ ^ i^^ rpj^g pressure and 
temperature are fixed at 1.0 atm and 300°C, respectively. 

E. Results 

1. Size Disparity Effects 

The effects of size disparity alone are presented in Fig. 1131 where we show the energy of mixing 
as a function of x = 0mi/0m for a blend (Mi = 10, M2 = 100) with en = 622 = £12 = —2.6 x 10^^^ 
J, so that £12 = 0. Thus, energetically, there is no preference. We take q = 14, and the pressure and 
temperature are fixed at 1.0 atm and 300°C, respectively. We show isobaric and DDIM results. 
While the energy of mixing is negative for the isobaric case, which does not correspond to isometric 
mixing, it is positive everywhere for DDIM, which does correspond to isometric mixing. This 
result should be contrasted with the Scatchard-Hildebrand conjecture © Q]) whose justification 
requires not only isometric mixing but also the London-Berthelot conjecture H10|) . We also show 
corresponding C12, which weakly changes with x. For AIi = M2 = 100, the energy and the volume 
of mixing vanish, which is not a surprising result as we have a symmetric blend (both components 
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identical in size and interaction). 

It is important to understand the significance of the difference in the behavior of AEm for the 
two processes in Fig. 1131 Prom Fig. [7| we find that the solubility parameter 5 at 1.0 atm and 25° C 
is about 9.844 and 9.905 (MPa)^'^ for M = 100, and 10 respectively. Thus, assuming the validity 
of the Scatchard-Hildebrand conjecture ©, we estimate AEm — 0.9 kPa at equal composition 
(x = 1/2). However, a correction for temperature difference needs to made, since the results in 
Fig. ^Jare for tc = 300° C. From the inset in Fig. [2 we observe that c almost decreases by a 
factor of 2/3, while their difference has increased at tc = 300°C relative to ic = 25°C. What one 
finds is that the corrected AEm is not far from the DDIM AEm in Fig. EJat x = 1/2, but has no 
relationship to the isobaric AEm, which not only is negative but also has a much larger magnitude. 

It should be noted that the pressures of the pure components in both calculations reported in 
Fig. ^Jis the same: 1.0 atm. Under this condition, the Scatchard-Hildebrand conjecture © cannot 
differentiate between different processes as 6i and 62 are unchanged. But this is most certainly not 
the case in Fig. 1131 For a symmetric blend, AEm = even in an exact theory. Since the symmetry 
requires 5i =62, we find that the London-Berthelot conjecture (|1()|) is satisfied. Thus, the violation 
of the Scatchard-Hildebrand conjecture we observe in Fig. ^Jis due to non-random mixing caused 
by size-disparity. 

2. Interaction Disparity Effects 

In Fig. 1141 we fix en = —2.6 x 10^^^ J, and plot the energy of mixing as a function of (—622) 
for a blend with Mi = M2 = 100. Thus, there is no size disparity but the energy disparity is 
present except when en = 622- We not only consider isobaric, and DDIM processes, but also 
the EDIM process. At en = e22, we have a symmetric blend; hence, AEm = and c^i = C22 
for all the processes. Correspondingly, we have C12 = cf'^or C22, and the correction I12 = for 
all the processes. Away from this point, AEm for the three processes are different, but remain 
non-negative. The energy disparity has produced much larger magnitudes of AEm than the size 
disparity alone; compare with Fig. 1131 This difference in the magnitudes of AEm is refiected in the 
magnitudes of C12, as shown in the figure. We again find that isobaric and DDIM processes are now 
quite different; compare the magnitudes of AEm in Figs. I13l and [ni The difference between DDIM 
and EDIM, though relatively small, is still present, again proving that isometric mixing alone is 
not sufficient to validate the Scatchard-Hildebrand conjecture. [We note that AEm can become 
negative (results not shown) if we add size disparity in addition to the interaction disparity.] The 
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FIG. 14: Energy of mixing, mutual cohesive energy density, and correction Z12 as a function of |e22| for a 
50-50 blend. We take Mi = Ah = 100, en = -2.6 x lO^^i j^ y^ == 1.6 x IQ-^s m^ and q = 14. The pressure 
and temperature are fixed at 1.0 atm and 3GG°C, respectively. 

non-negative AEm is due to the absence of size disparity, and is in accordance with the Scatchard- 
Hildebrand conjecture. Since AEm is the highest for DDIM, the corresponding C12 is the lowest. 
Similarly, the isobaric energy of mixing is usually the lowest, and the corresponding C12 usually 
the highest. We note that all ci2's continue to increase with |e22| • We observe that isobaric and 
EDIM ci2's are closer to each other, and both are higher than DDIM C12. In the inset, we also 
show li2, and we learn that it is not a small correction to the London-Berthelot conjecture H10|l 
for the isobaric case (non-isometric mixing). For isometric mixing, we have the usual behavior: a 
small correction /12. 

3. Variation with P 

We plot AEm, ci2) and /12 as a function of P in Fig. ^J for isobaric, DDIM, and EDIM 
processes. We note that all these quantities have a weak but monotonic dependence on P over 
the range considered. Again, the isobaric AEm remains the lowest, and consequently isobaric C12 
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FIG. 15: Energy of mixing, mutual cohesive energy density, and correction I12 as a function of pressure for 



a 50-50 blend. We take Mi = Ah = 100, en = -2.2 x lO^^^ J, 622 
and q = 14. The temperature is fixed at 300°C. 
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remains the highest over the range considered. As noted above, isobaric and EDIM ci2's are closer 
to each other, but different from DDIM C12 over the entire range in Fig. 1151 The monotonic 
behavior in P is correlated with a monotonic behavior in /12. We observe that I12 provides a small 
correction at 300° C over the pressure range considered here. The interesting observation is that 
isobaric /12 provides the biggest correction, while DDIM I12 the smallest correction. However, the 
EDIM li2 remains intermediate, just as EDIM C12 is. 

4. Variation with T 

We plot Ai^M and /12 as a function of T in Fig. ^Jfor isobaric, DDIM, EDIM processes along 
with the isochoric process. The corresponding Ci2's, all originating around 90 MPa at — 100°C 
as a function of T, are shown in Fig. 1171 All processes start from the same state at the lowest 
temperature (— 100°C) in the figure. There is a complicated dependence in AE'm on T over the 
range considered for some of the processes. Let us first consider the isochoric process in which all 
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FIG. 16: Energy of mixing, mutual cohesive energy density, and correction I12 as a function of temperature 
for a 50-50 blend. We take Mi = M2 = 100, en = -2.2 x lO^^i J, 632 = -2.6 x lO^^i J, i;o = 1.6 x 10"28 
m"^ and q = 14. The pressure is fixed at 1.0 atm. 

quantities show almost no dependence on T. The energy of mixing remains positive in accordance 
with the London-Berthelot conjecture (jlUj) . This is further confirmed by an almost constant C12, 
and an almost vanishing correction /12 in the inset. Both isobaric and EDIM processes give rise 
to negative AEm, thereby violating the London-Berthelot conjecture. The DDIM AEm shows 
a peak, which is about eight times in magnitude than its value at the lowest temperature, but 
remains positive throughout. The corresponding C12 shows a continuous decrease to zero with 
temperature for isobaric, DDIM, and EDIM processes. However, it is almost a constant for the 
isochoric process. The non-monotonic behavior in temperature of A£^m is correlated with a similar 
behavior in /12. This behavior is further studied in the next section. From Fig. 1161 we observe that 
/12 provides a small correction at 300° C. However, at much higher temperatures, it is no longer 
a small quantity, and depends strongly on the way the mixture is prepared, even if it remains 
isometric. In particular, isobaric I12 seems to provide the biggest correction. It is almost constant 
and insignificant for the isochoric and EDIM processes. 
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FIG. 17: Different mutual pressures as a function of temperature; system as described in Fig. E| 



F. Using Internal Pressure 



We now introduce a new quantity as another measure of the mutual cohesive pressure by using 
the internal pressure. To this end, we consider the internal pressure Pin in the RMA limit. We 
find from (0^ that 



RMA 



Pmvo > q[-eu4>^i/2 - e22(/>mi/2 - ei20mi<Am2], 



which can be expressed as 



RMA p 2/2 //P2 



PmVo ' CiiX (/'m/(/)ml + ^22(1 " x) (pml <Pm2 + 2ci2x(l - x), 

where we have used ()39() in the RMA limit. Now we take this equation as a guide to define a new 
mutual cohesive energy density, the mutual internal pressure P^ in the general case as 

P^ = [Pin - cX^x^cfi/ctCi - 4(1 " xfclt/<C2\/Ml - x). 



In the RMA hmit, Pj^ reduces to the RMA ci2 given in ^^. We show P^ in the inset in Figs. UHl 

)M 
IN 



[TCI andEl What we discover is that isochoric P^ is almost a constant as a function of temperature 
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around 120 MPa, and is larger than isochoric C12, which is around 90 MPa. Indeed, over most of 

IN 



the temperatures at the lower end, we find that PM > c\2- However, the main conclusion is that 



both C12 and P^ are monotonic decreasing or are almost constant, and behave identically except 
for their magnitudes. 



VI. NEW APPROACH USING SRS: SELF-INTERACTING REFERENCE STATE 

Unfortunately, the van Laar-Hildebrand cohesive energy cyi does not have the required property 
of vanishing with ei2, see (|39p . The unwanted behavior of cyi is due its definition in terms of the 
energy of mixing from which we need to subtract pure component (for which ei2 may be thought 
to be zero) cohesive energies c|^. The subtracted quantity is used to define C12, and if this definition 
has to have any physical significance, it should vanish in the hypothetical state, which we have 
earlier labelled SRS, in which eyi vanishes even though en and 622 are non-zero. The hypothetical 
state obviously violates the London condition ()11() . even though the real mixture does not. Let 
us demand the subtracted quantity to vanish for SRS. To appreciate this point, consider (|S7|) for 
SRS. It is clear that Ai^M continues to depend on the thermodynamic state of the mixture via ^a^ 
this quantity in general will not be equal to pure component quantity 0?^. With the use of H38|) in 
(|H7|) . we find that 

^E'^^Vq = eii[(/.ii - X(t)n,4>u/4>l,l] + 622 [fe " (1 " x)</.m022/'/'m2]> (47) 

which is usually going to depend on the process of mixing. On the other hand, from (|35|) . we 
observe that 

AEl^^'^'='{cf, + cl,M^2, (48) 

had ci2 = 0. In this case, AEm would depend only on pure component quantities cf^, and its 
behavior in a given process at fixed composition should be controlled by the behavior of cf^ in that 
process. This is not the case, as can be seen in Fig. I18| where we plot AEm for the hypothetical 
state SRS for a 50-50 mixture. We have ensured that the SRS state at the initial temperature 
in Fig. ^lis the same in isobaric and isochoric processes. But the pure components are slightly 
different. A new process is also considered in Fig. 1181 in which we set the volume Vsrs of the 
hypothetical SRS at a given temperature to be equal to the volume V of the real mixture (nonzero 
612) at 1.0 atm at that temperature. This process is labelled isobaric (I^rs = V) in the figure. 
The pure components are also at 1.0 atm for this process. Since the pure components for this 
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FIG. 18: A£;m for IRS as a function of temperature for a 50-50 blend. We take Afi = M2 = 100, 
en = -2.2 x lO^^i j^ ^^2 = -2.6 x lO^^i j ^^ ^ 1.6 x lO^^s m^ and q = 10. The pressure is fixed at 1.0 
atm. 

process are the same as in the isobaric calculation, (|48() requires AE'm for the two processes to be 
identical at all temperatures. This is evidently not the case. Consider Fig. 1111 From this, we see 
that c^ are almost constant with T for the isochoric case. Thus, according to (|15|). AE'm should be 
almost a constant, which is not the case. For the isobaric case, cj^ are monotonic decreasing with 
T, which will then make AE'm; according to H48|) . monotonic decreasing with T, which is also not 
the case. Hence, we conclude that the van Laar-Hildebrand C12 does not vanish for SRS. 

A similar unwanted feature is also present in the behavior of P^ introduced above for the same 
reason: it also does not vanish for SRS. 

It is disconcerting that the van Laar-Hildebrand 012(612 = 0) does not vanish for SRS, even 
though the mutual interaction energy ei2 = 0. This behavior is not hard to understand. The 
mixture energy is controlled by the process of mixing, since the mixture state varies with the 
process of mixing even at ei2=0. 
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SRS 



A. Mutual Energy of Interaction c^| 

To overcome this shortcoming, we introduce a new measure of the cohesive energy that has the 
desired property of vanishing with ei2- Let E denote the energy per unit volume of the mixture. 



We will follow 15|, and compare it with that of the hypothetical reference state SRS. Its energy 
per unit volume i^sRS differs from E due to the absence of the mutual interaction between the two 
components. Let Vsrs denote the volume of the SRS. In general, we have 

^SRS/^ = <Am/(Am,SRS- 

The difference 

E^^,^ =E- (ysRs/^)i?SRS 

represents the mutual energy of interaction per unit volume due to 1-2 contacts. From H20|) . we 
obtain 

-^int '^0 = eii[(/>ii - (^SRS/'^^)<?^11,SRS] + 622 [</'22 - (VsRs/^)</'22,SRs] + ei2(/>12, 

where the contact densities (pij without SRS are for the mixture state and with SRS are for the 
SRS state. These densities are evidently different in the two states but approach each other as 
ei2 -^ 0. The above excess energy should determine the mutual cohesive energy density, which 

„SR! 
"12 



we will denote by c^2 in the following to differentiate it with the van Laar-Hildebrand cohesive 



energy density ci2- It should be evident that Ef^^ vanishes as ei2 vanishes due to its definition. 

The absence of mutual interaction in SRS compared to the mixture causes SRS volume to 
expand relative to the mixture. This is shown in Fig. ^Jin which the void density in SRS is larger 
than that in the mixture. The relative volume of mixing at constant pressure (1.0 atm) for the 
mixture is shown in the upper inset, and that for SRS is shown in the lower inset. It is clear that 
the latter relative volume of mixing is always positive, indicating an effective repulsion between 
the two components. This should come as no surprise since the excess interaction ei2 for SRS from 

® 

^?P = -(eii + e22)/2>0, 

and has the value 2.4 x 10^^^ J in this case. Thus, it represents a much stronger mutual repulsion 
than the mutual repulsion due to ei2 = 0.1 x 10^^^ J in the mixture. The effect of adding the 
mutual interaction ei2 to SRS is to add mutual "attraction" that results in cohesion, and in the 
reduction of volume. Thus, the change in the volume can also be taken as a measure of cohesion, 
as we will discuss below. 
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FIG. 19: Void density and volume ol mixing as a function of temperature for a 50-50 blend. We take 
Ml = A/2 = 100, en = -2.2 x lO^^i j^ p^^ = _2.6 x lO^^i j^ ^^ ^ ^ g ^ IQ-^s ^3 ^^^ ^ ^ ^q rpj^g pressure 

is fixed at 1.0 atm. 



1. RMA Limit of E. 



(M) 



In the RMA limit, along with the monomer equality assumption ()4U|). which implies that 
1 SRS ) it is easily seen that 



j;^(M) RMA II j,f ,f n \ 



aP j,P 



AM) 



As a consequence, E^^^ /2x{l — x) reduces to the RMA value of (— C12) given in (|5U)) . This means 
that the quantity 0^2 defined via 



^SRS _ p(M) ,„ 



(49) 



has the required property that it not only reduces to the correct RMA value of the van Laar- 
Hildebrand C12, but it also vanishes with 612- [As usual, we assume the identification H41|) .] 
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FIG. 20: Mutual cohesive energy density, and correction I12 using the IRS as a function of temperature for 
a 50-50 blend. We take Mi = Afa = 100, en = -2.2 x IQ-^i J, 622 
and q = 10. The pressure is fixed at 1.0 atm. 
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2. New Mutual Cohesive Energy Density: c^l 



SR,S 



We now take H49() as the general definition of a more suitable quantity to play the role of cohesive 
energy density. This quantity is a true measure of the effect produced by mutual interaction energy, 
and which also vanishes with 612- Away from the RMA limit, C12 and c^2 '^^^ ^^o^ going to be 
the same. It is evident that C12 depends not only on T, P, or T, V, but also on the composition 
X = 4>mi/<Pm, Q, and the energies Cij. For isochoric calculations, we will ensure that the mixture 
and SRS have the same monomer density (j)^. In this case, Vsrs = V. For isobaric calculations, we 
will, as usual, ensure that they have the same pressure, so that the two volumes need not be the 
same. However, we will also consider Vsrs = ^ foi" isobaric calculations to see the effect of this on 
C12 ■ This process is what we have labelled isobaric (Vsrs = V) in the Fig. 1181 We show 0^2 ^o^ 
various processes in Fig. 1201 We see that it continues to increase monotonically for the isochoric 
case, and reaches an asymptotic value of about 50 MPa, while van Laar-Hildebrand C12 in Fig. 
Iltji is almost constant and about 90 MPa. An increase in solubility with temperature at constant 
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volume is captured by 0^2 but not by ci2- The behavior of the two quantities for the isobaric case 
is also profoundly different. While C12 monotonically decreases with temperature, cf^^ is most 
certainly not monotonic. It goes through a maximum around 400°C before continuing to decrease. 
This suggests that the solubility increases before decreasing. This should come as no surprise as 
we explain later in the following section. 

Note that while the definition of cf^^ is independent of any approximate theory, the definition 
of C12 depends on a form pSf) whose validity is of questionable origin as it depends on RST. It 
should be noted that the value of c^2 truly represents the energy change that occurs (in SRS) due 
to the presence of the mutual interaction between the two components. The changes brought about 
due to ei2 7^ should not be confused with the changes that occur when the two pure components 
are mixed ( ei2 7^ 0), because the latter involves not only the process of mixing when ei2 = (this 
gives rise to SRS ) but also involves the rearrangements of the two components due to the mutual 
interaction. It is the latter rearrangement due to mutual interaction that determines 0^2 • This 
is easily appreciated by studying the void density profiles for SRS and the mixture in Fig. 1191 
As usual, ei2 for the mixture is defined by pij) . The hypothetical process of mixing (under the 
condition ei2 = 0) also contributes to the energy of mixing, since the energy of this state is not 
the same as the sum of the two pure systems' energies, because SRS is affected by mixing the pure 
components. This is easily seen by computing the energy of mixing Ai^MjSRS for SRS for which 
ei2 = but not 611,622: 

VsrsAE'm.srs = VsRs-EsRS -Vi El - V^ E2 , 

where V^ ,Ef are the volume and energy density of the pure ith component. This part determines 
the energy of mixing due to pure mixing at ei2 = 0, and but it has nothing to do with the mutual 
interaction 612- Thus, this contribution should not be allowed to determine in part the mutual 
cohesive energy density. It is easy to see that 

AE'm = -Kjnt + (^srs/^)A£'m,srs, 

in which the second contribution is not a true measure of the mutual interaction. Nevertheless, it 
is considered a part of the van Laar-Hildebrand cohesive energy density. Using (|35() and (|49() in 
the above equation, we find that we can express c^2 i^ terms of the van Laar-Hildebrand cohesive 
energy densities as follows: 

c?P(ei2) = ci2(ei2) - ci2(ei2 = 0) + (1 - </'m/</'m,SRs)(cn + 42)/2- (50) 
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Here, cf2 (^12) is the SRS-based cf^^ introduced in (|1U)) . and 012(612) is energy-of-mixing based C12 
in (|c{5|) . In the special case (/)m = 0m,SRS; the last term vanishes, and the SRS-C12 is the difference 
of the van Laar-Hildebrand-ci2- In any case, because of the above relation, studying van Laar- 
Hildebrand-ci2 also allows us to learn about the SRS-ci2- Therefore, we have mostly investigate 
C12 in the present work. However, it should be realized that we also need van Laar-Hildebrand-ci2 
for SRS, data for which ate not available. 

The arguments similar to the presented above suggest that we can similarly use the difference 
^int— ^int,SRS) where Pint, SRS is the interaction pressure in SRS at the same volume V, as a measure 
of the cohesive pressure for the system. Alternatively, we can use SRS at the same pressure P as 
the mixture, whose volume Vsrs is different from V. We can use the negative of the difference 
V — Vsrs as a measure for cohesion. However, we will study these quantities elsewhere and not 
here. 

B. Cohesion and Volume 

We have argued above that the changes in the volume can also be taken as a measure of cohesion. 



We now fo 
elsewhere 



low_this line of thought. The volume of mixing is governed by two factors, as discussed 
: (i) the size disparity between the two components, and (ii) the interactions. To 



disentangle the two contributions, we can consider the difference 

AathV^ = V- Kth, 

where V^th is the volume of the athermal system (no interaction) at the same T, P. (In the 
athermal limit, AVM,ath = V^th ~^i ~ ^2 '^i^ be identically zero if the two polymer components 
have identical sizes (Mi = M2), and will be negative if they are different.) The difference Aath^ 
is governed only by the presence of interactions (en, 622, and 612). A positive A^th^ will imply an 
effective repulsion, and a negative AathV' will imply an effective attraction. Thus, AathV' can also 
be used as a measure of the cohesiveness of the mixture. 
It is easy to see that 

Vi,th/V = (kuj (t>m,a.th = (1 - </'o)/(l - </'0,ath), (51) 

in terms of total monomer density or void density. Using this, we can also calculate the difference 
Aath^^ = Aath^/^ per unit volume in our theory, and is shown in Fig. [^ in an isobaric process 
at 1.0 atm. We note that eoi = ^02 = £12 = (en = ei2 = 622 = 0) for the athermal state, 
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while £01 = 1.1 X 10-21 J, eo2 = 1-3 x IQ-^^ J, and £12 ^ 0.01 x lO'^^ J (en = -2.2 x lO'^^J, 
ei2 = —2.6 X 10^21 J, ei2 — —2.39 x 10^^^ J) for the mixture. Recall that we have decided to define 
ei2 for the mixture by (jllj) . The difference A^th^ is obviously controlled by all three energies eij, 
and not just by ei2. Since the mixture has attractive interactions (eij < 0), its volume V is smaller 
than Vath so that Aath'W < as seen in Fig. 1^ 

To determine the effect of only ei2, we use SRS in place of the athermal state and introduce 
the difference 

AsRSt" = 1- (^SRS/^), 

which is also shown in Fig. 1211 Here, the mixture is compared with the SRS state in which eoi = 
1.1x10-21 J, £02 = 1.3x10-21 J, butei2 = 2.4x10-21 J (en = -2.2 x 10^21 J, ei2 = -2.6xl0-2ij, 
ei2 = 0). Since the mixture has extra attractive interaction (ei2 < 0) than the SRS state (ei2 = 0), 
V is smaller than I^rS) so that the negative value of Asrs'W should not be a surprise. The minimum 
in AsRst) is due to the behavior of the void density, which is shown in Fig. ^J for SRS and the 
mixture. In the two insets, we also show the relative volume of mixing at constant pressure (1.0 
atm) for SRS (lower inset) and the mixture (upper inset). We notice that below around 500° C, 
the free volume in SRS changes faster than in the mixture, while above it, the converse is the case. 
This changeover in the rate of change causes the dip and the minimum in Asrs?^. Similarly, the 
dip in Aath"^ is caused by the changeover in the rate of change of free volume in the mixture and 
the athermal state. 

The minimum in Asrs^^, see Fig. 1211 also explains the intermediate maximum in the isobaric 
0^2 ill Fig. 1201 The behavior of c^2 shows that the mutual cohesiveness increases and then 
decreases with the temperature, and is merely a reflection of the way the void densities in Fig. 1191 
behave with the temperature. 

VII. MORE NUMERICAL RESULTS: 712 AND /12 

We observe from (|39|) that the quantity 

ci2Vo/i-qei24'l,/2) 

approaches 1 in the RMA limit. As we have seen earlier (compare Figs. l^andfTU]). the quantity q/2 
in the RMA limit should be properly replaced by {q/2 — u) for finite q. The new factor incorporates 
the correction due to end groups. Therefore, we introduce a new dimensionless quantity 

712 = -ci2VQ/{q/2 - v)ei24>'i,, 
51 
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FIG. 21: Volume difference Aath^ and Asirsw for a 50-50 blend as a function of temperature. We have 
Ml = M2 = 100, en = -2.2 x lO^^i j^ ^^^^ ^ _2.6 x lO^^i j y^ = 1.6 x lO^^s ^3 g^^j ^ ^ ^0. Tfie pressure 
is fixed at 1.0 atm. 

and investigate how close it is to 1. Any deviation is due to nonrandoni mixing and will provide a 
clue to its importance. In Fig. 1221 we plot 712 as a function of temperature for various processes. 
The coordination number is taken to be g = 14, and the initial void density is 0o = 0.00031333. 
This system is identical to the one studied in Fig. 1161 For the isochoric process, 712 is almost a 
constant and very close to (but below) 1. It decreases gradually for EDIM and becomes almost 
0.92 at about tc = 2000°C. The situation is very different for isobaric and DDIM processes where 
712 shows enhancement at a temperature near the boiling of one of the pure components. To be 
convinced that the peak in 712 in Fig. |521is indeed due to boiling, we plot it for g = 14 in Fig. 
1231 We see a discontinuity in 712 for all cases except for isochoric case, as is evident from the inset 
(in which isobaric results are not shown). The discontinuity is due to the boiling transition which 
occurs at different temperatures for the different cases. These discontinuities become rounded as 
in Fig. 1221 when we are close to boiling transitions. 

We now turn our attention to the binary correction I12, which measures the deviation from the 
London-Berthelot conjecture (|10|) . The existence of this deviation in the case when we choose ei2 
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FIG. 22: The scaled quantity 712 for various processes. This quantity is very close to 1 except at high 
temperatures where one of the components seems to be near its boiling. For the system considered, we have 
Ml = M2 = 100, en = -2.2 x lO^^i j^ ^^^ = _2.6 x lO^^^ J wo = 1-6 x IQ-^s ^3 ^^^^ ^ ^ i^_ j'^^ pressure 
is fixed at 1.0 atm. 



to satisfy the London conjecture ()11() will suggest that it is caused by thermodynamics which makes 
cohesive energies different from their respective van der Walls energies e^. Moreover, the behavior 
of /12 also reflects partially the behavior of the pure components, since the definition of /12 utilizes 
the pure component cohesive energies, see ()12() . even if C12 is defined without any association with 
them as is the case with 0^2 , whose definition is independent of pure component quantities c|^. 
Consider the inset in Fig. 1141 which shows I12 associated with ci2- We observe that it is almost 
constant for the two cases (DDIM and EDIM) for which mixing is isometric. However, it has a 
strong variation for the isobaric case where mixing is not isometric. Thus, we see the first hint 
of a dependence on volume of mixing in the behavior of li2- Such a dependence on the volume of 
mixing is also seen in the inset showing /12 (corresponding to C12) in Fig. 1161 We are struck by the 
presence of maxima and minima for the two isochoric cases. They appear to be almost at same 
temperatures in both cases; notice the location of minima and maxima in the isobaric /12 and the 
volume of mixing. 
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FIG. 23: The scaled quantity 712 for various processes for a 50 — 50 mixture. This quantity is very close 
to 1 except at high temperatures where the sharp discontinuity is due to one of its components undergoing 
boiling. The system is Mi = 100, M2 = 10, en = -2.6 x lO^^^ J= 622, vq = 1.6 x lO'^^ m^ and q = 14. 
The pressure is fixed at 1.0 atm. 

To further illustrate the relationship between isobaric /12 and the volume of mixing Aum, we 
plot them as a function of the temperature in Figs. |^ and |531 for q = 6 and q = 14, respectively, 
for several values of the pressure P. The strong correlation between their maxima and minima 
are clearly evident. These figures illustrate that the binary correction I12 for the isobaric process 
closely follows how the volume of mixing Aum behaves. A negative At;M is a consequence of the 
effective attraction between the mixing components, which corresponds to a larger value of their 
mutual cohesive energy density C12; this is reflected in a corresponding negative value of li2- As 
the pressure increases, the correction decreases and the system gets closer to satisfying the the 
London-Berthelot conjecture H1U() . Thus, we can conclude that the binary correction is negligible, 
but not zero, as long as the volume of mixing is very small or even zero. 
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FIG. 24: Binary deviation Z12 and the relative volume of mixing in the inset for a 50 — 50 mixture. For the 
system considered, we have Mi = M2 = 100, en = -2.2 x IQ-^i j^ 622 = -2.6 x lO^^i j i,^ ^ 1.6 x IQ-^s 
m^ and q = 6. 

VIII. DISCUSSION AND CONCLUSIONS 

We have carried out a comprehensive investigation of the classical concept of cohesiveness as 
applied to a thermodynamic system composed of linear polymers. In pure components, we are only 
dealing with one kind of interaction: microscopic self-interactions between the monomers. There- 
fore, we are interested in relating this microscopic self-interaction to a thermodynamic quantity 
in order to estimate the strength of the microscopic self-interaction. In binary mixtures (polymer 
solutions or blends), there are two self interactions, and one mutual interaction. The microscopic 
self interactions in the mixture do not depend on the state of the mixture. Therefore, they are 
evidently the same as in pure components. Thus, it is the microscopic mutual interaction that is 
new, and needs to be extracted by physical measurements. However, what one measures exper- 
imentally is not the microscopic interaction strengths, but macroscopic interaction strengths due 
to thermodynamic modifications. 

Traditionally, cohesive energies have been used as indicators of the macroscopic interaction 
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FIG. 25: Binary deviation I12 and the relative volume of mixing in the inset for a 50 — 50 mixture. The 
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strengths. The definition of c for pure components, given in ((2)), does not depend on any particular 
approximation or any theory. It is a general definition. Thus, it can be measured directly if the 
interaction energy iSint can be measured. It can also be calculated by the use of any theory. Only 
the obtained value will depend on the nature of the theory. 

On the other hand, the mutual cohesive energy density C12 is defined by a particular form of the 
energy of mixing. This gives rise to two serious limitations of the definition. The first deficiency 
is caused by the use of the form (|H5|) . whose validity is questionable. Thus, the extracted value 
of C12 is based on this questionable form. The second deficiency is that this value also depends 
on pure component properties c^, since the mutual energy is measured with respect to the pure 
components. However, we have argued that the energy of the mixture will be different from the 
sum of the pure component energies even if ei2 is absent. Thus, C12 does not necessarily measure 
mutual interactions. One of our aims was to see if some thermodynamic quantity can be identified 
that could play the role of a true C12 that would vanish with 612- This is the quantity c^2 that we 
have identified in this work. There is no such problem for c^, since it not only vanishes with e, as 
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we have seen before, but it also does not suffer from any approximation. 

We follow our recent work [l5| in the approach that we take here. We confine ourselves to a 
lattice model, since the classical approach taken by van Laar __!] and Hildebrand [3] is also based 
on lattice models in that their theory can only be justified consistently by a lattice model. We then 
introduce two different reference states to be used for pure components and mixtures, respectively. 
For pure components, we introduce a hypothetical reference state in which this self-interaction 
is absent. Since there is only one interaction, which is absent in the reference state, this state is 
nothing but the athermal state of the pure component. With the use of this state, we introduce the 
interaction energy <Sint defined via (0) . Because of the subtraction in (Q) , Si^t depends directly on 
the strength of the self-interaction, the only interaction present in a pure component that we wish 
to estimate. Thus, it is not a surprise that i^int is an appropriate thermodynamic quantity that 
can be used to estimate the strength of the self-interaction. Consequently, we use this quantity 
to properly define the pure component cohesive energy density c at any temperature. We have 
discussed how this definition resembles the conventional definition of c in the literature. For the 
mixture, we introduce the self-interacting reference state (SRS), which allows us to define c^2 
that vanishes with ei2. In contrast, the van Laar-Hildebrand C12 does not share this property. 

As noted earlier, our approach allows us to calculate c^ and C12 even for macromolecules like 
polymers, which is of central interest in this work or for polar solvents, which we do not consider 
here but hope to consider in a separate publication. The cohesive energy density c^ is a macro- 
scopic, i.e. a thermodynamic quantity characterizing microscopic interparticle interactions in a 
pure component. In general, being a macroscopic quantity, c is a function of the lattice coordi- 
nation number q, the degree of polymerization M, and the interaction energy e, in addition to the 
thermodynamic state variables T, and P or V. We have investigated all these dependences in our 
study here. The same is also true for C12. 

Based on the same philosophy, we need to introduce another hypothetical reference state called 
the self-interacting reference state (SRS) for a blend, which differs from a real blend in that the 
mutual interactions are absent, but self-interactions are the same. The difference £^^ introduced 
in Q allows us to estimate the strength of the microscopic mutual interaction between the two 
components of a blend; we denote this estimate by cf^^ to distinguish it from the customary 
quantity C12 originally due to van Laar, and Hildebrand and coworkers. However, £j^^ is not what 
is usually used to define the mutual cohesive energy density C12. Rather, one considers the energy 
of mixing AEm- We have compared the two approaches in this work. 

The conventional van Laar-Hildebrand approach to solubility is based on the use of the regular 
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solution theory. Thus, several of its consequences suffer from the limitations of the regular solution 
theory (RST). One of the most severe limitations, as discussed in the Introduction, is that the 
theory treats a solution as if it is incompressible. This is far from the truth for a real system. 
Hence, one of the aims of this work is to investigate the finite compressibility effects, i.e., the 
so-called "equation-of-state" effects. As a real solution most definitely does not obey the regular 
solution theory, the experimental results usually will not conform to the predictions of the theory. 
Thus, we have revisited the solubility ideas within the framework of a new theory developed in our 
group. This theory goes beyond the regular solution theory and also incorporates equation-of-state 
effects. We have already found that this theory is able to explain various experimental observations, 
which the regular solution theory (Flory-Huggins theory for polymers) cannot explain. Therefore, 
we believe that the predictions of this theory are closer to real observations than those of the 
regular solution theory. 

A. Pure Components 

We have first studied the cohesive energy density c of a pure component. Since it is determined 
by the energy density of vaporization, it is a quantity independent of the temperature. It is also 
clear from the discussion of the van der Waals fluid that c , defined by , is independent of the 
temperature, since the parameter a in (|.S0)) is considered T- independent. This particular aspect of 
c^ ensures the equality of Pin = —-Pint and {88 / dV)T\ see (|S5|) . For a pure component, Hildebrand 
J2] has argued that the solubility of a given solute in different solvents is determined by relative 
magnitudes of internal pressures, at least for nonpolar fiuids. Thus, we have also investigated the 
internal pressure. However, the equality ()33() is not valid in general as we have shown above. Thus, 
the internal pressure, while a reliable and alternative measure of the cohesion of the system in its 
own right, cannot be equated with {d£/dV)T- Their equality is shown to hold only in RMA. If we 
allow a in pOf) to have a temperature-dependence, the equality H33|) will no longer remain valid. 
Moreover, it can be easily checked that under this situation. Pin and (? for the van der Waals fiuid 
will no longer be the same. 

We have investigated the pure component c using our recursive theory. We find that it is 
very different depending on whether we consider an isochoric process, in which case it is almost 
constant with T, or an isobaric process, in which case it gradually decreases to zero at very high 
T. In general, we flnd that 

Cy > Cp. 
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We have found that isochoric c and Pin are almost constants as a function of T. It provides a 
strong argument in support of their usefulness as a suitable candidate for the cohesive pressure. But 
c^ and Pin do not remain almost constant in every process, as the isobaric results in Fig. EJ clearly 
establish. The same figure also shows in its inset that the isobaric c^ exhibits a discontinuity due 
to boiling, as expected. Unfortunately, most of the experiments are done under isobaric conditions; 
hence the use of isochoric cohesive quantities may not be useful, and even misleading and care has 
to be exercised. We also see that c changes a lot over the temperature range and replacing it by 
its value at the boiling point may be not very useful in all cases. 

The cohesive energy also depends on the molecular weight, and usually decreases with M as 
expected; see Fig. |7| This is the situation for q = 14, and should be contrasted with the behavior 
shown in the inset in Fig. [3 where we show c^ for M = 10, and 100, but for q = 10. All other 
parameters are the same as in Fig. |7| What we observe is that at low temperatures, c^ for M = 10 
lies above the Cp for M = 100, while the situation is reversed at high temperatures. This crossover 
is due to the boiling transition that the smaller M pure component must undergo at about 600° C. 
Thus, Cp decreases with increasing M only far below the boiling temperatures. The dependence 
on the lattice coordination number q is also not surprising: it increases with q but requires end 
group correction in that c^ is proportional to (g/2 — v); see Fig. El This means that the solubility 
function 6 also depends on the process, and its value at the boiling temperature of the system will 
be dramatically different in the two processes. We have also found, see Fig. El that c^ is not truly 
a quadratic function of the monomer density, one of the predictions of the regular solution theory. 

B. Mixtures (Solutions or Blends) 

We now turn to the mutual cohesive energy density ci2- Here, the situation is muddled since 
the definition of ci2 is based on a form (|,'-{5j) of the energy of mixing AEm, whose validity is 
questionable beyond RST. This should be contrasted with the definition of c, defined by (jSJ, which 
is independent of any particular theory. Therefore, c^ can be calculated in any theory without any 
further approximation except those inherent in the theory. It can also be measured directly by 
measuring iSint, and does not require any further approximation to extract it. On the other hand, 
the calculation of ci2 in any theory requires calculating AEyi in that theory; thus, this calculation 
is based on the approximations inherent in the theory. However, one must still extract ci2 from 
the calculated AEm by expressing AE'm in the form (|35jl . As discussed above, this form is justified 
only in RST and ci2 in 1)351) at this level of the approximation is indeed a direct measure of the 
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mutual interaction energy ei2- Whether ci2 defined by the form (|35|) still has a direct dependence 
on ei2 is one of the questions we have investigated here by introducing SRS. For c\2 to be a direct 
measure of the mutual interaction energy ei2, we have to ensure that it vanish in SRS. What we 
have shown is that c\2 does not vanish in SRS as seen in Fig. 1181 A new measure of the mutual 
cohesive energy density c^^ is introduced in (|49|) . which has the desired property. 

The reference state SRS behaves very different from the mixture or its athermal analog, as 
clearly seen from the Figs. 1181191 and 1211 In particular, SRS has strong repulsive interactions 
between the two species. This results in the SRS volume to be much larger than that of the 
mixture. Therefore, it is possible that the two components may undergo phase separation in SRS. 
Its use to define (^^ ., therefore, should only be limited to the case where SRS is a single phase 
state so as to compare with the mixture. We do not report any result when SRS is not a single 
phase. 

We have found that c^^ has the correct behavior that it first rises and then decreases with 
T for a compressible blend. As explained above, this is a reflection of the way the void density 
behaves with the temperature, as shown in Fig. 1211 The rise and fall of cohesion is not apparent in 
the temperature-dependence of c\2- It has also been shown that c^^^ can be expressed in terms of 
Ci2. Therefore, we have basically explored the behavior of c\2 more than that of (^^ ■ As expected, 
and as shown in Fig. I12L c\2 increases with q. However, it is not just simply a linear function of q^ 
as seen in Fig. 1121 It is also clear from the behavior in this figure that cyi changes its curvature 
with T. Thus, cyi is not linear with the inverse temperature /? over a wide range of temperatures. 

We have been specifically interested in the contribution due to volume of mixing. For this reason, 
we have considered three different kinds of isometric mixing (zero volume of mixing) , two of which 
we have called EDIM and DDIM. These are introduced in Sect. V. The third isometric mixing 
is studied as part of isochoric processes, which we have also considered. The fourth process is an 
isobaric process in which mixing is at constant pressure. We have also compared cyi with a related 
internal pressure quantity P^ , which is another measure of the mutual cohesive energy density or 
pressure. For a particular process, the two quantities behave similarly, though their magnitudes 
are different in that P^ > c\2\ see Fig. 1171 We have found that isobaric and EDIM quantities are 
somewhat similar, and both are different from the DDIM quantities. All three quantities have a 
strong temperature dependence, but the isochoric quantity is very different in that c\2 or P^ for 
the latter remains almost a constant with T. 

We have paid special attention to the violation of the Scatchard-Hildebrand conjecture (jSJ even 
when the microscopic interactions obey the London conjecture ()lip. We find that even under 
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isometric mixing (EDIM), the energy of mixing can become negative, as seen in Fig. 1161 There 
is another possible source of violation of the Scatchard-Hildebrand conjecture seen in Fig. 1131 
Here, the pure components in the isobaric and DDIM processes remain the same, so that the pure 
component solubilities S are the same. Therefore, from ^ we would conclude that the energy of 
mixing should be zero, regardless of the process. What we find is that the energy of mixing is not 
only not zero, it is also different in the two processes. This violation is due to non-random mixing 
caused by size and/or interaction strength disparity. 

We have also investigated the behavior of the binary correction I12. We find that it need not 
be small. In particular, we find that it becomes very large in isobaric and DDIM processes as 
seen in Fig. 1161 As the pressure increases, /12 decreases in magnitude and the deviation from the 
London-Berthelot conjecture (jlUj) becomes smaller. This effect shows that the deviation from the 
London-Berthelot conjecture is due to the presence of compressibility. Thus, the behavior of the 
free volume determines that of li2- In the isobaric case, we have discovered a very interesting fact. 
The behavior of /12 mimics the behavior of the volume of mixing as seen clearly in Figs. 1241 and l251 
This observation requires further investigation to see if there is something unusual about isobaric 
processes or it is jut due to non-isometric mixing in the presence of free volume.. 

In summary, we have found that it is not only the non-isometric mixing by itself that controls 
the behavior of the cohesive quantities, but also the non-random contributions. 
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